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SYNOPSIS 

DYNATIC EESPOISE 0"^ COlPOSITE PLATES Y/ITH 
CUT-OUTS ANL HOLES 

A Thesis Submitted 

In Partial Fulfilment of the Eequirements 
For the Degree of 
3X)CT0E OF PHILOSOPHY 

by 

A. RAjArani 
to the 

Department of Mechanical Engineering 
Indian Institute of Technology, Kanpur 
March, 1978 

The study of the dynamic behaviour of structural elements ma.de 
of composite materials has attracted considerable attention in recent 
years. In spite of the many recent investigations in the area of 
vibration of composite plates, not much attention has been focussed 
on the dynamic response of composite plates with discontinuities. 

The thesis first reviews, in Chapter 1, the literature on the study of 
dynamic response of composite plates. The principal objective of the 
present work is to develop a general computation technique based on 
Lagrange's equations of motion and to study the effect of disconti- 
nuities on the free vibration and blast response of isotropic and 
composite plates. The work can be broadly divided into three sections; 

1) theoretical study on the free vibration characteristics' of a 
composite plate with discontinuities, 

2) experimental verification of these characteristics on a 
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3) investigation into the effect of blast loading on the response 

of composite plates. 

The free vibration characteristics of a composite plate with central 
cut-outs and circular holes is studied respectively in Chapters 2 and 3. 

The effect of discontinuity is assumed to be siiailar to a displacement- 
dependent loading. Lagrange's equations of motion have been employed 
and these lead to an infinite system of frequency equations in the case 
of free vibrations. A suitable size of the frequency equations is chosen, 
striking a compromise between the accuracy and the computation time 
involved. The convergence of the results has been studied. Square and 
rectangular plates are caisidcred as examples. Two edge conditions, 
namely simply-supported and clamped-clamped, are taken for detailed 
invest igation . 

The analysis of composite plates with cut-outs is simpler. The 
integrals necessary in the development of frequency equations have been 
derived in closed form for both simple-support and clamped-clamped edg’O 
conditions. For clamped-clamped plates, these integrals are given in an 
Appendix. However, in the analysis of plates with holes a numerical 
integration scheme has been employed to evaluate these integrals. The 
results are first compared with available results on isotropic solid 
plates and it is found that the results obtained are sufficiently 
accurate. A good agreement is also found with the available results on 
isotropic plates with cut-outs and holes. The cause of discrepancy, if 
any, has been explained. Numerical results have been computed for 
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square^ simply -supported and clamped-clampod composite plates for 
different fibre-orientations. Ratios of principal Young's moduli of 1 
(balanced bidirectionally reinforced glass-epoxy), 3(unidirectionally 
reinforced ^ass-opoxy), 10(unidirectionally reinforced boron-epoxy) 
and 40 (unidircctionally reinforced graphite-epoxy) have been used in 
the computation. The behaviour of frequencies in the case of balanced 
bidirectional plates with cut-outs is found very similar to isotropic 
plates and hence no results have been computed for these plates with 
holes. The effect of various paramoters such as the ratio of the cut-out 
dimension (or hole dimension) to plate length, fib re -orientation angle 
and modulus ratio is censidered on the natural frequencies and mode shapes 
of the composite plate. Some interesting observations have been made on 
the behaviour of natural frequencies and mode shapes. 

The second part of the work is concerned with the experimental 
verification of theoretical results obtained earlier. Chapter 4 describes 
the experimental set-up and procedure. Only square plates with centrally 
located circular holes are considered lin the present experimenta.l study, 

A test fixture was designed and fabricated to conduct experiments on the 
free vibration studies of composite plates with holes. The simple-support 
and cl amped-clamped boundary conditions were simulated by first conducting 
experiments an isotropic solid plates and verifying the results with 
available theoretical results. A fairly good agreement showed that these 
edge conditions have been approximated. The effect of clamping arrangement 
hos also been taken into account. 
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The unidirectionally reinforced glass-epoxy composite plates were 
made by a hand lay-up technique. The ijHterial was characterized to 
determine its static and dynamic properties. The dynamic properties Imve 
been used in the computation of numerical results. Experiments were 
conducted on simply- supported and clamp ed-clarnped square plates Y/ith 
central circular holes. The frequencies and mode shapes have been deter- 
mined as the hole size was enlarged in steps in the aluminium plate as 
well as in the unidirectionally reinforced fibre-glass epoxy plates with 
fibre -orient at ions .0^3 15^? 30^ and 45^* The influence of the holes on 
the natural frequencies and mode shapes have been shown. There is a good 
agrecimcnt between the theoretical results and experimental results. A 
comparison of theoretical results for the hole and the cut-out is also 
made and it is found that small cut-outs and holes influence the frequencies 
in a similar way. The mode shapes are found to bo influenced due to the 
presence of cut-outs and holes. 

Chapter 5 deals with the last phase of the work - the transient 
response of composite plates v\^ith and without holes. The modal approach 
has been used in the computation of niumorical results. The results have 
been obtained for isotropic and unidirectionally reinforced E-glass-epoxy 
platos and are given for tv\^o cases s l) V7hen there is no hole in the plate 
and 2 ) when there is a centrally located hole in the plate. Results have 
been obtained for both simply-supported and clamped-clamped boundary 
conditions but arc given only for clamped-clamped plates. In order to 
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verify the theoretical results, experiments have been conducted on isotro- 
pic and inai directional E-glass-epoxy rectan.gular plates using a shock tube 
as the loading device. The expeiimental peak dynamic strains have been 
normalized with respect to the pressure and corupared with theoretical peak 
strains. The strain- time histoiy is also shOYm for a particular gauge 
location in the unidirectional reinforced glass-epoxy plate. Einally, a 
comparison, of djmamic amplification factors (theoretical) ha.s been made 
between an isotropic and a unidirectional glass-epcxj"- plate. 

In Chapter 6, some general conclusions of the present work are 
sumriBiized and possible extensions of the present work are discussed. 



CHAPTEE 1 


IITRODUCTION 


1 . 1 G-eneral 

The potential use of fibre-reinforced composite materials for 
various design applications in aerospace, marine and terresterial structures 
is well known. An improvement in the performance of these structures 
requires a study of the dynamic response characteristics of the structural 
elements, the most common being plates and shells, made of composite 
materials. Structural and architectural requirements may necessitate the 
presence of discontinuities either geometrical or material in these 
structural elements. Geometrical discontinuities such as cut-outs or holes 
may lighten the structuie , alter the resonant frequency and are needed for 
venting purposes. The present investigation is concerned with the dynamic 
response of uni directionally reinforced composite plates with the geometri- 
cal discontinuities. 

In the last two decades considerable amount of work has been done on 
the dynamic response of composite plates. The work done on the dynamic 
behaviour of these plates can be placed into four groups. 

In the first, the plate is studied using ah approximate theory based 
on Kirchoff hypotheses specifying the behaviour of the plate through its 
thickness with shear deformation neglected. The theories based on Kirchoff 
hypotheses have been developed by Pister and Dong CO. Eeissner and 
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Stavsky CO Dong, Pister and Taylor , Stavsky Ql-I] summarized 

in [52 . 

In the case of high degree of orthotropy in the plate, an extension 
of Mindlin’s theory for homogeneous isotropic plates CO has been proposed 
by Yang, Norris and Stavsl 5 y CtI and by Whilney and Pagano CbH. In these 
theories linear displacements are again assumed across the plate thickness. 
Transverse shear deformation, however, is not neglected. Whitney DI 
obtained closed form solution for bending deflection, flexural vibration 
frequencies and buckling loads of simply-supported rectangular plates of ; 
special construction and showed the significance of shear deformation on i 
the response of hi^ly anisotropic laminated plates. Sun and V/hitney DoD 
also showed the effect o f heterogeneous shear deformation over the thickness 
of the plate on the dynamical behaviour of general laminates. These 
theories are classified in the second group. 

The third category may be characterized as approximate elasticity 
solutions in which attempts were made to bridge the gap between 

the approximate plate theories and the general theory of elasticity formula- 
tion by smoothing and averaging procedures. These approaches lead to 
'Effective Modulus' and 'Effective Stiffness' theories which are again 
approximate* 

In the final group, vibrations of an (infinite) laminated plate are 
based on the complete three-dimensional theory of elasticity where the 
solutions for all layers are combined in accordance with the interface 
traction and displacement conditions to give a frequency equation. 
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However, for all but the simplest case of two layer con struct ion 
the algebra involved in both the generation and solution of frequency 
equations become intractable. 


1.2 Review of Previous \¥ork 

A fairly exhaustive review of the literature on the mechanics 
and dynamics of structures made of composite materials can be found in 
three seivey papers by Bert and Prancis Gap and Bert 09,2^1. 

Huffington and Hoppmann [21 Hiavestigated the free transverse 
vibration of thin rectangular orthotropic and isotropic plates for a variety 
of boundary conditions all of which have the characteristics that two 
opposite boundaries were simply-supported. Levi's method Q£] has bean 
used to obtain frequency equations and mode shapes. HearrnOn applied 

Rayleigh^ s method to derive closed foimulas for the vibration of specially 
orthotropic plates using beam characteristic functions appropriate to the 
boundary conditions. Dickinson (54^ extended a sine series solution 
developed by Dill and Pister 05]] to the vibration of orthotropic plates. 
Plates with two opposite edges simply -supported, clamped along all four 
edges 9 and clamped along two opposite edges and free along the other edges 
were taken as examples. Swarup and Sundararajan natural 

frequencies and buckling loads jof generally ortho tropic rectangular plates 
with different combinations of simply ^supported, clamped and free edges usjng 

Rayleigh-Ritz procedure. 
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Ashton and Tifaddoups [27^] presented a Rayleigh-Eitz solution for the 
buckling and free vibration of plane anisotropic rectangular plates with 
various boundary conditions. They obtained natural frequencies and 
mode shapes for free-free and contilever plates and these compared well 
with their experimental results. Whitney and leissa (^28529^ obtained 
closed form solutions for the case of arbitrarily laminated orthotropic 
plates with simple support boundary conditions. Their results showed a 
strong effect of b ending -st re ching coupling in lowering the frequencies. 

An analytical and experimental study of transverse vibrations of 
laminated orthotropic plates was carried out by Hikami []30^ involving 
combinations of simply -supported and free edge conditions. He considered 
1) all edges simply -supported, 2) three edges simply-supported , fourth free, 

3) two opposite edges simply-supported, other two free and 4) two adjacent 
edges simply -supported, other two free. The solution of the vibration 
problem for the first three edge conditions was obtained by the method of 
liOvy ( 2 ^ and for the fourth-edge condition was obtained by the Rayleigh-Ritz 
method using product of beam characteristic functions. A good agreement 
was obtained between the theoretical and experimental results for symmetric 
modes. 

Ashton and Anderson DO studied the natural .frequencies and mode 
shapes of clamped-clamped boron-epoxy plabes using Ritz method. They 
obtained ^cellent agreement between analytical and experimental results. 

Bert and Mayberiy [^2^ presented an analysis based on Rayleigh-Ritz method 
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to determine the natural frequencies of vibration of laminated cross— ply 

rectangular plates. They gave results for fully clanped boundary conditions 
and compared them with experimental results also. Ashton also studied 

the natural frequencies and characteristic mode shapes of free-free plates 
and illustrated the effect of 1) the orientation of the constituent plies, 

2) the stacking sequence, and 3) the degree of orthotropy. 

Mohan and Kingsbury [3^ made an analytical study of the natural 
frequencies and mode shapes of anisotropic thin plates by use of Galerkin's 
method. They investigated the change in natural frequencies and mode 
shapes as the material elastic axes of the plate were rotated with respect 
to the geometrical edges. 

Y/hitney [35 ,362] showed that use of energy methods in conjunction 
with classical beam characteristic function yields very slow convergence 
of the vibration frequencies for hi^ly anisotropic materials when the 
natural boundary conditions are not satisfied. He obtained improved 
results by using a classical double litourier Series L3i1 analysis in 
which both geometric and natural boundaiy conditions are satisfied. 

Decapua and Sun [383 obtained frequencies and mode shapes of a cl^s of 
orthotropic plates with either clamped or simply -supported boundaries 
using Eayl eigh-Eitz method. Maurizi and Laura [393 

frequencies of vibration of clamped rectangular plates by ■using a simple 
polynomial expression and applying Galerkin method. 
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Claiy and Cooper [40[[]used a finite element analysis to study the 
vibration characteristics o f cantilever wing like flat trapezoidal 
aluminium plates reinforced with boron-epoxy (placed symmetrically about 
the middle surface). The analytical results were in qualitative agreement 
with experimental trends for nodal patterns and frequency variations 
with filament orientation. Thornton and Clary correlation 

study of fimt e— element modelling for vibrations of composite material 
panels. They obtained excellent agreement of their results with those 
given by Ashton and Anderson [3^Zl • Rosettos and Tong [3C] hybrid 

stress finite element method to study the effect o f fibre-orientation 
angle and orthotropy parameter (ratio of Principal Youngs s moduli) on the 
vibration and flutter characteristics of cantilevered anisotropic plates, 

Jones presented an exact solution and numerical results for 

simply -supported plates that are laminated unsymme t ri ca lly about their 
middle surface. Elishakoff gave frequencies of clamped plates in a 

form analogous to the corresponding frequencies of a simply-supported 
plate* In the case of a simply -supported plate^ the v#ave numbers aie equal 
to mtr /a and ntr /a respectively and m and n are positive integers which 
determine the mode shape. Eor a clamped plate the wave numbers were 
presented in the fDrm pit /a and qtr /a where the pair of real quantities 
p and q could be found from the solution of two supplementaiy eigenvalue 
problems. Recently lin and King jj-SZl extended the work of Whitney and 
leissa m to the more general case of two opposite edges simply aipported 
and two with arbitrary edge conditions. 
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Thus a study of the vibrations of composite plates has been made 
earlier using double Fourier series method, Rayleigh-Rit z method, Galerkin 
method and the finite element method. Ro work seems to have been done on ! 
the vibrations of composite plates mth cut-outs or holes. Earlier work 
done include vibrations of isotropic plates and shells with cut-outs and/or 
holes. Kumai found three natural frequencies of vibrations of square 

plates with circular holes analytically as well experimentally. He computed 
the natural frequencies for a given mode by means of superposition of 
solutions of circular ring plates satisfying the boundary conditions at 
several points along the edges of the square plate. His approximation 
allowed small residual deflections at some portion of the outer boundary of 
the square plate. Takahashi analyzed the problem of rectangular 

plates with circular holes and all ends clamped using Rayleigh -Ritz method. 
Joga Rao and Pickett Qs^] also used the Rayleigh-Ritz method to find the 
fundamental frequency of vibrations of square plates with circular holes 
and square cut-outs. They modified the deflection functions for the mode 
shapes of solid plates to account fsr the appropriate singularity due to the 
presence of the circular hole. Brogan et.al. Z49J did a systeniatic 
analytical and experimental investigation of the dynamic behaviour of a 
cylinder with a .rectangular cut-out. They approximated the potential and 
kinetic energies of a thin shell by two-dimensional finite difference 
i^thods. They obtained excellent agreement between experimental and 
theoretical results. Anderson, Irons and Zienkieuicz computed 

natural frequencies of vibrations of clamped square plates with central 
circular holes using finite element method. They used triangular elements 
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and the mesh size was finer near the hole boundaiy. Results were given 
for the first and fourth mode frequencies only. Kristiansen and Soedel (3lZI 
found fjndamental frequencies of clamped, square plates with various types 
of cut-outs using Rayleigh Method. Paramsivam C523 described a method 
of determining the effects of openings on the natural frequencies with 
different types of boundaiy conditions using a grid frame work model. 
Basdekas and Chi natural frequencies of a simply-supported 

isotropic plate with square cut-outs assuming the effect of cut-out to be 
equivalent to a displacement-dependent loading. In the methodology used 
by them, a knowledge of the normal modes is necessaiy. They extended this 
approach to investigate the dynamic response of oddly reinforced shells I3C] 
and the non-linear forced vibrations of a non-prismatic beam [^553 • Aksu 
and Ali f 56 | predicted the <^namic characteristics of orthotropic plates 
by a finite difference unequal interval formulation and extended this 
formulation to determine the ^namic characteristics of rectangular plates 
with cut-outs [571. Further studies have been made on the vibrations of 
cylinders with cut-outs recently [58, 59, 6 o3 . 

The study of the transient response of structures is quite complex. 
Earlier studies were made on the elastic and plastic response of structural 
elements subjected to impact and shock loading. Elastic response studies 
have been made by Baker ct.al. []6l3» Cheng and Beneveniste [62 , 633 » 
Banerjee QO , Crocker [653, Crocker and Hudson 11663, and Rajamani and 
Sundararajan [^T3 elastic beams and plates. A good review of the work 
done on the dynamic plastic behaviour of plates and shells can be found 
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^ -A- few studies have also been made on the dynamic response 

of composite plates. These include the works done by Moon Ql^ ? Chow QCI? 

Sun and ¥/hitney [^73 9 74? 753 ? Sieu and Bert Dsn , Sim and Chattopadhyay Ddl 
and Kubo and Nelson [j^. 

1.3 Present Tfork 

The principal obiective of the present work is to develop a general 
— based on Lagrange's equations of motion and to study 
the effect of geometrical discontinuities on the free vibration and blast 

isotropic and composite plates with cut-outs and holes. The work 
can be broadly divided into three categories: 

1) theoretical study on the free vibration characteristics of a 
composite plate with discontinuities, 

2) experimental verification of these characteristics on a unidircc- 
tionally reinforced glass-epoxy plate with holes, 

3) investigation into the blast response of composite plates with and 
without holes. 

The free vibration characteristics of a composite plate with central 
cut-outs and circular holes is studied respectively in Chapters 2 and 3. The 

ps assumed to be equivalent to a displacement 
dependent loading. Lagrange's equations of motion have been employed and 
these lead to an infinite system of frequency equations in the case of 
free vibrations. A suitable size pf the frequency equations is chosen- 
strikiiig a compromiSG between the accuracy and the computation time involved. 
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Sq[uare plates with two edge conditions, simply -supported and clanped-clanped, 
are taken for detailed investi^tion. 


The analysis of composite plates with cut-outs is simpler* The 
integrals necessary in the development of frequency equations have been 
derived in closed form for both simple-support and clamped-clamped 
boundaiy conditions. Tbr clamped-clamped plates, these integrals are given 
in Appaidix A. In the analysis of plates with holes, the integrals cannot 
be evaluated in closed form and therefore, a numerical integration scheme 
■has been employed to evaluate these integrals. The results are first 
coupared with available results and it is found that the results obtained 
are sufficiently accurate. A good agreement is also found with the 
available results on isotropic plates with cut-outs and holes. The cause 
0 f di screpancy, if any, has been explained. Numerical resiolts have been 
computed for square, simply -supported and clamped-clamped composite plates 
for different fibre-orientation an^es. Ratios of principal Yoxmg's moduli 
of 1 (balanced bidirectionally reinforced glass-epoxy), 3 (imidirectionally 
reinforced glass-epoxy), 10 (unidirectionally reinforced boron-epoxy) and 
40 (unidirectionally reinforced graphite -epoxy) have been used in the 


computation. The behaviour of frequencies in the case of balanced bidirec- 
tional plates is found very similar to isotropic plates and hence no results 
have been computed for th^e plates with holes. The effect of various 
parameters such as the ratio of the cut-out dimension (or hole dimension) to 
plate length, fibre-orientation angle and modulus ratio is considered on the 
natural frequencies and mode shapes of the composite plate. Some interesting 
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observations have been emde on the behaviour of natural frequencies and ' 
mo de shap es . 

The second part of the work is concerned with the experimental 

verification of theoretical results obtained earlier. Chapter 4 describes 

the experimental set-up and procedure. Only square plates with centra,lly 

located circular holes are considered in the present experimental study, 

A test fixture was designed and fabricated to conduct experiments on the 

free vibration studies of composite plates with holes. The simple-support 

and clamped -clamped boundary conditions were simulated by first conducting 

eaperiments an isotropic solid plates and verifying the results with 
A 

available theoretical results. A fairly good agreement showed that these 
edge conditions have been approximabed. The effect of clamping arrangement 
has also been taken into account . 

The -uni direct ion ally reinforced glass-epoxy composite plates were 
made by a hand lay-up technique. The material was characterized to 
determine its static and dynamic properties. The dynamic properties have 
been used in the computation of numerical results. Experiments were 
conducted on simply-supported and clamped-clamped square plates with 
central circular holes. The frequencies and mode shapes have been deter- 
mined as the hole size was enlarged in steps in the aluminium plate as well 
as in the uni directionally reinforced fibre-glass epoxy plates with fibre- 
orientations 0°, 15°; 30° and 45°, The influence of the holes on the 
natural frequencies and mode shapes have been shown. There is a good 
agreement between the theoretical results and exp erimenta.1 results. 
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A comparison of theoretical results for the hole and the cut-out is also 
made and it is foiond that small cut-outs and holes influence the frequencies 
in a similar vray. The mode shapes are found to be influenced due to the 
presence of cut-outs and holes. 

Chapter 5 deals with the last phase of the work - the transient 
response of composite plates with and without holes. The modal approach 
has been used in the computation of numerical results. The results have 
been obtained fbr isotropic and xml directionally reinforced E-glass-epoxy 
plates and are given for two cases; l) when there is no hole in the plate 
and 2) when there is a centrally located hole in the plate. Results have f 
been obtained for both s imply -suppoiled and claraped-clamped boxmdary j "■ 
conditions but are given only for clamped -clamped plates. In order to 
verify the theoretical results, experiments have been conducted on isotropic 
and unidirectional E-glass-epoxy rectangular plates using a shock tube 
as the loading device. The experimental peak dynamic strains have, been 
normalized with respect to the pressure and compared with theoretical peak 
strains. The strain-time history is also shown for a particular gauge 
location in the xmidirectional reinforced glass-epoxy plate. Finally, a 
comparison of dynamic amplification factors (theoretical) has been made 
between an isotropic and a unidirectional glass-epoxy plate. 

In Chapter 6, some general conclusions of the present work are 
summarized and possible extensions of the present work are discussed. 



CHAPTEE 2 


VIERATIOUS OP COMPOSIIDE PIATES WITH CUT-OUTS 

2.1 Introduction 

One of the problems in structural dynamics is the vibrations of 
plates or shells with discontinuities. Several approximate methods, 
such as Eayleigh-Ritz finite difference [49, gi'id frame work 

analogy Qssjhave been previously used to study the dynamic behavio\jx 
of plates or shells with cut-outs. In the present chapter, natural 
frequencies of generally orthotropic plates with cut-outs are --determined 
by using Lagorange's equations of motion. The potential and kinetic 
energies based on Kirchoff's plate theory are used to dete.piiine the 
motion of the plate. The effect of the cut-out is expressed as a dis- 
placement-dependent loading similar to the treatment of Basdekas and 
Chi DU. The present approach, however, is more general and the mode 1 
shapes of the"referQace plate" need not be eigenfunctions as in 
reference [^53^ • In the methodology used by Basdekas and Chi a 
knowledge of the normal modes is necessary. The analysis eii 5 )loys an 
assumed space-mode approach with multimode expansions by taking charac- 
teristic beam functions and finally an infinite system of second order 
differential equations is obtained in generalized co-ordinates.. A 
suitable size of the frequency equations is chosen based on the accuracy 
and convergence of results. The mode of convergence has been discussed. 
Results are given for square, simply -supported and clamped-plamped 
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composite plates with square cut-outs for different composite mateidals. 

A comparison is made with available results on isotropic plates. Certain 

* 

interesting observations have been made and discussed in this <^pter. 

2,2 Equations of Motion 

The Lagrange's equations of motion in the presence of nonconserva- 
tive forces have the form 


_d^ 

dt 




3T 


3V 


== P 


r = 1,2,, 


( 2 * 1 ) 


where T is the kinetic energy, V is the potential energy, q^ are the 
time -dependent generalized displacement co-ordinates and P^ are time- 
dependent generalized forces. 

In case of plates with cut-outs, the generalized forces Pj,(t) will 
include the effect of cut-outs. 


The potential energy for a solid uniform plate, shown in Pighre 
2;1a due to bending is 


V 


^ gP "^0 


+ 2T)^-n- w , w + D w + 4D(i;c w ^ , 
,xx 1‘2 ,xx ,yy 22 ,yy 66 ,xy 


+ 4D^>' w w + 4D„> w ... w ■ “IdA 
1>6' ,xx ,xy 26 ,yy ,xy -i 


(2-;2) 


where w is the noriiial deflection of the- plate and a comma denotes pa^ial 
differentiation, 

* A part of Chaplfer 2 is reported in |^79,8o3* 
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"i3 's can be written as CaO 

^ J ~\- 1 ^ 

for each layer are given below in terms of the invariants [3 1]®’“-'^ 
fibre -orientation angle for the layer with respect to the geometrical 
axes of the plate. 

+ Ug cos 20+ cos 4 0 
Q22 = - Ug cos 20 + Uj cos 49 

S2 "" ^4 " % 

%6 " ^5 ~^3 

— 1 

Q^g = “ "2 ^2 ~ ^3 

_ 1 

Q 25 = “ "S’ ^2 ^ ® + Uj sin 4 6 (2.3) 

in which, 

Ui = (5Q^^ 30^2 + 2Q^2 ^ 

^2 = (^U - *^22^/^ 

U3 = (0,1 + Q22 - 2Qi2 - 4 Qgg )/8 

^4 ” ^ '*’ ^*^12 ~ 

U5 = (Qn 022 - 2 Qi 2 + 4Qgg)/8 
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The kinetic energy of the solid plate is given by. 


a b 

T = / J ph TO dj;‘i 


(2.5) 


where p is the mass density 


and h is the thickness of the plate. 


The effect of size, shape and location of the cut-out can be consi- 


dered as a di placement-dependent external loading and, therefore, it is 


expedient to express w as 


E S 

w(x,y,t) = I I q (t) 4 U,y) 


( 2 . 6 ) 


m=1 r!;=1 


where q (t) is the generalized co-ord-inate in the m mode^ 


4 (xov) are admissible fimctions in the mn mode, 

mn^ ^ - 


Substituting for in the expressions for T and Y in equation (2^1) 


and simplifying 3 
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+ <!>.. <i> ) + (<l> <l>. + 4> 6. )ZI<iA) q, 
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A ^ , 

3 ~ 1 j 2 j**.S 


where is the number of cut-outs, 
c 

Defining and P„. ., the above equation can 

® mn’ ran ’ mn’ mn P13 

be written as, 
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The terms with negative sign in equation (2.8) can be thought of as 
an external loading on the plate and these terms are present because of 
the cut-outs in the plate. The equation of motion can be written as. 


!, 1 V 

iit=1 n=1 
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(2.14) 


5t)r free vibrations, equation (2,14) reduces to 

in + i i - “m) V = 
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= 0 


i — 1,2,..., R 


( 2 . 15 ) 
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Equation (2.15) defines an eigenvalue problem and its solution 
gives the natural frequencies and corresponding eigenvectors. Ihe eigen- 
value problem has been solved for two bovindaiy conditions: l) simply- 
supported edge conditions, 2) clamped-clamped edge conditions. 

2.2,1 Simply -Supported Plates 

The simply -supported edge boundary conditions are, 


X = Oja ; 

11 

o 

II 

-hi ",xx 
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12 ,yy 

- 2D w 

16 ,xy 
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Ii 
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22 ,yy 

- 2D w 

26 ,xy 

0 


(2.16) 

It is not possible to give an exact solution for the eigenvalue 
problem. However, a good approximate solution can be obtained by ^suming 
that the function (j) ^ (x,y) can be given as a product of beam char act eristi* 
functions in x and y directions. 

For a simply -supported plate, 

w(x,y,t) = I I S sin ^ (2.17) 

m=1 n=1 

<[) (x,y) = sin sin (2.18) 

The function <|>jjj^(x,y) in equation (2.18) satisfies the geometric 
boundaiy conditions but not the natural boxmdaiy conditions of equation (2,16) 
except in the case of an isotropic plate. 

ii ii ii iri 

Substituting * in the expressiors ibr M , 3 , C and « 

° ^ mn’ mn’ mn mn 
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In the case of an Isotropic plate, the equation of motion (2.15) 
reduces to a simple form which is also given by Basdekas and chi tZ53 ^ as 

*^ij “ig ‘^ij "" [fkij ~ J I ^mn ^ “ I I Om 

(2.25) 

This simplification is possible only when 4* (x>y) are eigen- 

mn 

funotionb in the mn mode, 

2,2,2 Clamped-Clamped Plates 

The boundary conditions for a clamped-clamped plate are given as 
x=0, a; w = 0 , w = 0 

(2.26) 

ysOob: w=0 , w = 0 

’ ,7 

A good approximation to the mode shape of the plate can be obtained 
by assuming the function «(i^(x,y) to be a product of beam characteristics 
functions in x and y directions , 

(x,y) = X^(x) Y^(y) 

where , 

X (x) = cosh - cos A (sinh sin ) (2,2?) 

a am' a a 

Y^(y) is also given by a similar expression. The calculated values of 
and Ajjj are given in Table 2.1 . 

» 

ConsMering only one cut-out, the expressions for and 

mn' '^mn' mn 

equations (2.9 - 2.12) can be obtained by defining the following 

integrals , 
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where a comma denotes partial differentiation. 

= ^a/=* > = V® ! 

and Xg are the distances of the cut-out edges from the origin and are 
shown in figure 2 . 1 b . The integrals in equations ( 2 . 28 ) and ( 2.2 9 ) 

are given in Appendix A. 

The expressions for and are 


\n ^ ^ i) d) 


( 2 . 30 ) 
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( 2 . 31 ) 
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Table 2,1 


Parameters for a clamped- clamped beam 


m or n 


Prequencgr parameter a or a. 


m 


n 


Modal coefficient A or A 

m n 


1 

4.730040744862703 

0.9825022145762379 

2 

7.853204624095836 

1 .000777311907269 

,3 

10.99560783800167 

0.9999664501254087 

4 

14.13716549125746 

1.000001449897656 

5 

17.27875965739948 

0.9999999373443833 

6 

20.42035224562606 

1 .000000002707595 

7 

23.56194490204045 

0.9999999998829942 

8 

26.70353755550818 

1 .00000000005056 

9 

29.84513020910325 

0.9999999999997815 

10 

32.98672286269281 

1.000000000000009 

11 

36.12831551628262 

0.9999999999999996 
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Unlike in the ease of an isotropic simply -supported plate, the 
equation of motion (2,15) can not be reduced to a simple form for a 
clamped-clamped plate as are not eigjnfunctions , Therefore, one 

cannot start with a differential equation of motion as given in 
references 135, 543 


2,3 Computation technique 

The method of solution in deriving equations of motion assumed the 
presence of discontinuities in plates as displa cement -depdn dent external 
loading* An effective use of this method involves evaluation of numerous 
definite inte^^s necessary in the expressions for M^ , C^ and 

ot^« The coii 5 )utational technique employed product of beam characteristic 
functions for the admissible functions as tables of integrals [^8 2 3 
available and the definite integrals in equation (2.29) can be evaluated 
in closed form without much difficul'^. Once these integrals are evaluated, 
they can be stored in the computer memory and are accessible whenever needed* 


The natural fireguencies and corresponding normalized eigenvectors have been 
determined using the power method* This is an iterative method with 
orthogonalization. The method has been found to give very accurate 3 ?esults 


for certain known cases. 


2.4 Results and Discussion 

The analysis presented in Section 2.2 has been applied to sqpaare, 
simply— Supports d and elanped -clamped plates with centrally located square 
cut-outs. Natural frequencies and eorrespondii^ eigenvectors were obtained 
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employing power method with orthogonalization for different fibre 
orientations in steps of 15^* Four groups of symmetric ^ antisymmetric 
and a combination of symmetric and antisymmetric modes are considered. 
Ratios of principal Young ^ s moduli of and 40 have been used in 

the computation. In the case of uni directionally reinforced composite 
plates results are given for fibre orientation angles 0^ and 45^ only 
in terms of non-dimensional freq^uency^ o) , 

The material properties of the plates for obtaining the numerical 
results are given in Table 2.2. 

2.4.1 Mode of Convergence 


The Lagrange *s equations of motion, for free vibrations, lead to 
an infinite system of second order differential equations in generalized 
co-ordinates. It is essential to choose a proper size of the system of 
equations and to check the accuracy and convergence of frequency values 
so obtained. For this purpose, two cases are considered. For an isotropic 
simply -Supported plate, a comparison is made between the frequency values 
Obtained by -taking 16, 25 and 36 equations for a combination of symmetric 
modes. It is found that the convergence of the frequency values is very 
good. The frequency values and the number of iterations required to get 
them for the same size of the system are given in Table 2.3. The 
frequency values do not differ upto seven significant digits for 
c/a = 0,0. Twenty five equations are taken for all the cases censidered. 
Once the size of the system is fixed, a compailson in the frequency values 
is made again by taking ^irmBtric, aatisymmetric and a combination of 
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Table 2.2 


Material properties used in the computation. 




Property 


Mat erial 


G ^ 2/®2 

^2 

Isotropic 

1 

0*385 

0.30 

Balanced 

Bidirectional 

glass-epo3ty 

1 

0.200 

0.10 

Uni dire ctional 
glassy-epoxy 

3 

0*5 

0.25 

Unidirectional 

boron-epoxy 

10 

0.333 

0.30 

Unidirectional 
graphite -epoxy 

40 

0.500 

0.25 
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The numbers in the paranthesis show the number of iterations required in the 
eigenvalue evaluation. 

The numbers in the paranthesis correspond to frequency values for the new mode. 
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symmetric and antisymmetric modes. Two values of c/a (O.O and 0.3) are 
considered. It is observed that the frequenqy values for a particular ^ 
mode do not differ by taking any combination of modes for a solid isotropic 
plate (c/a = O.O) . This is found to be t lue for an orthotropic solid plate 
and also for clamped boundary conditions. The above observation, however, 
is not found true for a plate with a cut-out. When a comparison, between 
the frequency values obtained by taking all modes and other combination 
of modes is made, the values are found to differ (see Table 2.5). The 
maximu.m error is found in the second mode for an isotropic simply -supported 
plate (the results for other cases are not given) and is approximately 
7.5 percent. The error, however, is considerably less in symmetric modes. 

(for example, first, fourth and seventh in Table 2.3). A relative compari- 
son in .the frequency values led to the decision to compute natural frequencies 
by taking 25 equations and different combination of modes for all the cases 
considered. The different combination of modes is preferred as it is 
easier to keep track of the modes .and the number of iterations required l^o^ 
compute the natural frequencies in the solution of eigenvalue problem is i 
reduced. 

It is worthwhile to mention here that the a_ccuracy of the results 
depends to an appreciable extent on the capability of the employed admissible 
functions to approximate the true shape of the deflected surface to the 
required degree. An increase in the number of terms tends to improve the 
accuracy of the results in ^neral, but the resulting increase in the size 
of the matrices leads to increased computation costs and the possible 
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round-off errors may limit the accuracy of the results. 

2.4.2 Simply -Supported Plates 

a) Effect of cut-out on the natural frequency 

Table 2.4 gives a comparison of the non-dimensional frequencies for 
cut-out parameters equal to 0,0 and 0*5 in the case of an isotropic plate. 

The fundamental frequency obtained by Paramsivam lower than the 

exact value whereas the present metho d yiel is exact values. This may be 
due to the fact that a coarse mesh size has been taken in Ref. DO , Even 
the value of the fundamental frequency obtained by Aksu and Ali |^56 ^ 

solid isotropic plate using finite difference scheme with optimum unequal 
intervals is lower than the exact value. A comparison of first four modes 
is also made with the available results of Paramsivam Basdekas [^53 ]] 

for a plate with c/a equal to 0,5. The results obtained by Basdekas D3 n 
for the fundamental frequency is approximate as the influence of other 
modes has not been considered. When the influence of other modes is taken 
into consideration an extremely good agreement is obtained between the 
frequency values for two eases considered. The results are given in 
Table 2.5 for two values of c/a. This agreement is due to the reason 
that the equations of motion in the present method and in the approach 
used by Basdekas and Chi [^5321 same. 

The variation of frequencies with cut-out parameter, e/a, for an 
isotropic plate is shown in Figure 2,2. The frequencies have been 
subscripted (11, 12, 21, etc.) on the basis of largest amplitude coefficient 
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Table 2.4 

Comparison of non-dimensional frequency for a simply -supported 
isotropic plate for c/a = 0.0 and c/a = 0,5- 


c/a 

0.0 


0 

.5 


Mode 

Source 

1 

1 

2 

3 

4 

Param^ivam 

19,630 

25.450 

42.110 

70.500 

78.650 

C520 






Basdekas 

19.739 

25.082 

- 

- 

- 

C53l 






Present 

19.739 

23.627 

41.634 

72.134 

80.229 

Exact 

19.739 

- 

-- 

- 

- 

Aksu and 






Ali 

19.719 


— 

- 

- 

11563 







Table 2.5 


Comparison of sjrmmetric and all modes in terms of frequency ratio 
(= 4a)/“^^) for a simply -sup ported isotropic plate, c/a =0.2 

and c/a = 0.6 


Mode 

Symmetric Modes 

All Mode's 


c/a = 0.2 

c/a = 0.6 

c/a = 0.2 

c/a = 0.6 

1 

3 .8605 

8.3811 

4.0551 

8.5042 


(3.86f 

(8.35) 

(4.05) 

(8.50) 

2 

94.197 

57.629 

25.051 

23.214 


(94.2) 

(57.5) 

- 

- 

3 

111.59 

231.51 

61 .059 

64.631 


(111) 

(231) 

(61.0) 

(64.5) 

4 

346 .89 

1153.0 - 

94.551 

79.135 


(347) 

(835) 

(94.5) 

(78.5) 

5 

— 

— 

114.56 

(114) 

159.49 

6 



167.82 

246.44 

(246) 


a. The results in parentheses are from Reference {^552] 
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in the corresponding normalized eigenvectors. A fairly good agreement is 
obtained for the fundamental frequency with the results obtained in DO 
and in D3l. There is a sli^t decrease in the fundamental frequency 
for small cut-outs but the frequency increases as the cut-out size is 
increased. The behaviour of frequencies ( “ ^21^^ *^22* ^23 ^32^’ 

0 )^^ is found to be similar. They first decrease with the inescease in cut- 
out size but increase as the cut-out size is increased. The frequency, 

, increases with the increase in c/a but at larger c/a ratios, it 
decreases abruptly. Two interesting features have been observed^ (l) the 
existence. of a new mode due to the presence of a cut-out and (2) the 
cross-over of frequencies. It is seen that the fourth mode is identical 
with the fifth mode for a solid isotropic plate. The presence of a cut-out 
gives rise to the fourth mode. This is observed for frequencies correspond- 
ing to both syinmetric and antisymmetric modes *^24’ ^26 

not for frequencies corresponding to a combination of symmetric and anti- / 

, . r- — — etc.). The new mode, for example the fourth 

symmetric modes “ 23 ’ ^14 

mode, has largest equal amplitude coefficients A^j and with opposite 
signs whereas the fifth mode has largest equal amplitude coefficients 
and A_^ with same sign. The intro (iict ion of a small discontinuity seems 
to perturb the system for v/hich an exact solution is known and gives rise 
to a new mode. 

The other interesting feature observed is the cross-over of modes. 
Modes corresponding to and and these corresponding to and ix>^^ 
interchange* The mode shapes at the cross-over and nearby points have been 
obtained and shown later for an orthotropic plate. It is observed that the 
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mode shape for a particular mode at the cross-over and nearby points does 
not charge. Moreover; at the cross-over frequencies the plate has the 
mode shape corresponding to either of the inodes which cross over. 


In Figures 2.3 and 2.4; the variation of m with c/a is shown for 
0=0° and 0 = 45° respectively; for balanced bidirectionally reinforced • 
plates. It is seen that decreases with c/a for 0 = 0° and then 
increases for larger cut-out sizes. However; Ihere is no decrease in the 
value of the fundamental frequency with increasing c/a values for 0 = 45°. 
A similarity between the balanced-bidirectional plate and the isotropic 
plate is observed. In the case of balanced bidirectional plates also the 
modes corresponding to and are the same for a solid plate and the 
presence of a cut-out gives rise to a new mode corresponding to frequency 

^ 3 - 


It is seen in Figure 2.3 that the mode corresponding to inter- 
changes twice with that corresponding to 9 = 45°; the mo do 

corresponding to interchanges first with that corresponding to 

and then with that corresponding to There is a tendency for the 

modes corresponding to andw^^ to interchange and for the modes 


corresponding to 1^22 “15 


to separate from each other (after interchanging 


once) at higher c/a values as 0 is increased. The frequencies 


0) 


12 


and 


coincide for all fibre orientation angles in the balanced bidirectional 
case. 


The variation of m with c/a is shown in Figures 2.5 and 2,6 for 
imi directionally reinforced glass fibre-epoxy plates for 0 = 0° and 
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9 = 45° » respectively. For 6 = 0 °, decreases very little with 

increasing c/a values and then increases as the cut-out size is increased. 
There is no decrease in with c/a for 6 = 45° just as in the case of 
balanced bidirectional plates. It is seen that as the fibre orientation 
angle is increased, curves for ti )^2 closer together and when 

0 = 450 , they coincide. This is also true, for ixi and for small c/a.- 

ratios. In this case also some of the modes interchange: the mode 
corresponding to interchanges with that corresponding to when 

0 = 0 °; mode corresponding to interchanges with those corresponding 

to a).j^ and when G = 45°. However, modes corresponding to and 

interchange with each other once and twice for fibre orientations 0 ° and 45° 
respectively. To study the behaviour of mode shapes at the cross-over and 
near by points, these have been obtained for and cOgg shown in 
Figure 2.7 for fibre-orientation angle equal to 0°. It is observed that ( 
the mode shape for a particular mode at the cross-over and nearby points I 

does not change the plate can have either of the mode shapes at the cross- 

' — / ■ ' " ■ ■ ■"' 

over point. As in the case of balanced bidirectional plates, when 0 is 
increased, there is a tendency for the modes corresponding to <^22 ^ 3 -] 

to interchange and for the modes corresponding to 0^22 “-53 

from each other after interchanging once at higher c/a values. 

In Figures 2.8 and 2.9, the variation of o) with c/a is shown for 

unidirectional boron-epoxy plates. It is seen (5'igure 2,8) that all 

frequencies except wT decrease with increase in c/a values and at higher 

•3 

cut-outs they increase for 6 = 0 °; the rate of decrease being maxiiiium for 



39 


^21* However^ for 0 = 45^9 there is no decrease (figure 2#9)* 

Unlike isotropic 9 balanced bidirectional and unidirectional glass-epo:35r ; 
platos, there is no interchange of modes corresponding to (^^2 ^‘13 

for 0 = but in the case of 0 2= 45^9 these modes separate from each 
other after interchanging once at higher c/a ratio as in the previous 
cases considered; 

for unidirectionally reinforced graphite-eposy plates ^ figures 2*10 
and 2J1 show the variation of O) with c/a for 6 = o'^and 6 = 45^. When 
6 = 9 0)^^ and decrease with increasing cut-out sizes and the rate 

of decrease is more as compared to isotropic, bidirectional glass-epo^y, 
unidirectional glass-epoxy and unidirectional boron-epoxy plates. However, 
for 0 =r 45^9 (0^^ increases with increase in cut-out size. It is interesting 
t o not e, in figure 2 .10 that even the modes corresponding to and 
interchange for small (c/a = 0.3) cut-out sizes, but they exhibit normal 
behaviour for e/a greater than about 0.5* Moreover, for 0 = 0^ , the mode 
corresponding to ^oes not interchange with that corresponding to 

whereas they are found to interchange for balanced bidirectional and 
unidirectional glass-epoxy plates. 

At this stage it is necessaiy to comment on Ihe behaviour of 
frequencies. It is generally foimd that frequencies corresponding to 
different modes behave in different ways. However, some similarity is 
observed in this behaviour. The fundamental frequency for an orthotropic 
plate ( 6z= 0®) is found to decrease for medium cut-out parameters 
(c/a 2r 0.2 - 0.5) unlike in the case of isotropic plates. This indicates 




c/a 


f?ig,2.2 Ip vs c/a t^^ sirn ply-supported is 



w/v Don/Pha 
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that cut-outs make orthotropic plates less stiff for medium cut-out sizes 
in all modulus ratios considerod. Frequencies “^2’ ^21’ ^22^ ^23 ^14 

are fcund to decrease with increase in cut-out size for all fibre orientation 
angles and modulus ratios. They increase vi?ith larger cut-out ratios. The 
behaviour of frequency is, however, unpredictable. 


An attempt has been made to study the behaviour of frequencies. For 
this, the deflection at the centre is found in the case of a unidirectionally 
reinferced glass-epoxy, simply -supported plate which gives an idea about the 
variation in stiffness (and hence the fundamental frequency) with the cut- 
out ratio. In Figure 2.12 which shows the variation of the deflection v?ith 
cut-out ratio for a unidirectionally reinforced glass-epoxy plate and for | 

I 

fibre orientation an^es equal to 0 ° , 15°, 50° and 45°, the stiffness is ! 

found to decrease initially with c/a| the rate of decrease being maximum j 

! 

for 0=0°. As c/a is increased (c/a >0.4), stiffhess starts increasing i 
for all the fibre orientation angles. This explains the decrease in the 
fundamental frequency at medium cut-outs and the increase for larger cut-outs. 
The explaaation of the behaviour of higher modes is, however, dif fi cult . 


b) The effect of fibre-orientation an^e on the natural frequency 


The effect of fib re -orientation an^e on the frequencies and 
has been shown for different values of cut-out . ratios. In Figures 2.13 
to 2.16, the variation of with 6 for different c/a values is shown for 

balanced bidirectional, unidirectional glass-epoxy, unidirectional boron- 
epoxy and luii directional graphite -epoxy plates. It is seen that in all the 
cases increases with 6 for all cut-out parametfeJs,* aI the' cutug^ size is 
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increased, tlie rate of increase of increases, For balanced bidirectional 
plates (Figure 2.13), corresponding to e/a = 0.2 is lower than that 
corresponding to c/a = 0.0 at 6 = 0° and the curves for corresponding 
to these two values of c/a intersect as 8 is increased. In Figures 2.14a 
and 2.141, a ■unidirectional glass-epos^'" plate, is shown to vary 

with 6 for different cut-out parameters, the variation of with 0 for 
smaller cut-outs being shown in Figure 2.14a. It is observed in Figure 
2,l4a that is minimum for c/a = 0.2 at 0 = 0°. Ihe frequency, however, 
increases with 8 . In Figure 2.14b, it is found that the frequency curves 
for c/a = 0.0 and 0.2 intersect just as in the case of balanced bidirectional 
plates. For unidirectional boron-epoxy and unidirectional graphite -epoxy 
plates (Figures 2.15 and 2.16 respectively), the frequency curves for 
c/a = 0.0 and c/a = 0.2 do not intersect and the frequency values become 
almost equal at 8 = 45°. However, in Figures 2.15 and 2.16 it is seen 
that the curve for c/a = 0.0 intersects with curves for c/a = 0.3, 0.4, 0.5 
and 0.6 in case of unidirectional boron-epoxy plates whereas it is found to 
intersect with these for c/a = 0.3, 0.4, 0.5, 0.6 and 0.7 in case of uni- 
directional graphite -epoxy plates. Thus, there is a decrease in the 
fundamental frequency with increase in c/a in the case of simply-supported 
plates with higher modulus ratios for small values of ^ . However, there 
is noj^end in this decrease. 

The variation of different c/a values is shown in 

Figures 2.17 to 2.20 where it is foimd to increase with 0 for all cut-out 
parameters similar to It is also observed that for a given value of 0 , 















r different Guffouf pdrametiefs 
y-supported plate; R = J O 
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the frequency decreases with increase in c/a ve.lues whereas this 

behaviour was observed in case of only for siaaller values of . The 
behaviour of other frequencies, for example 0)22 etc.^have not been 

showi. Mohan and Kingsbury p~ 3 ^also observed the increase in<»J^^, 
and 0 ^ 2 ^ with fibre-orientation angle for a square, simply -supported 
boron-epoxy plate with no cut-outs. This trend is also found here for 
all the modulus ratios considered. 

c) The effect of modulus ratio on the natural frequency 

In order to understand the effect of modulus ratio on the natural 
frequency of plates with cut-outs, the behaviour of for different 
modulus ratios is shown with c/a in Figure 2.21 in the case of fibre 
orientation angles 0° and 45°. It is observed that for isotropic plates, 
decreases for small cut-out ratios whereas for an orthotropic plate 
(0 = 0°), the decrease in occurs upto larger values of c/a. The rate 
of initial decrease is found to increase with increase in modulus ratios. 

For 0 = 0°, the curves fox- the fimdamental frequency in the case of 
E^/E 2 = 3 and E^/E 2 = 10 are found to intersect a.t a value of c/a = 0.638 

approximately. This suggests that the increase in orthotropy may not be 
useful for plates with higher cut-out ratios when the fibre-orientation 
angle is 0° as stiffness is considerably reduced. It is also seen that I 
the fLindamental firequency for a balanced bidirectional orthotropic glass- 
epo3qy plate is always smaller than that for an isotropi c plate for all 
cut-out parameters. The variation of with c/a for different modulus 
ratios in the case of 0 = 45° is sjxown by dotted lines. The fundamental 
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frequency is found to increase for all modulus ratios with c/a* Also, the 
fuin da. mental frequency for a balanced bidirectional glass-epo:}^ plate is 
greater than tha.t for an isotropic plate* Also, the frequency curves for 
different modulus ratios do not intersect. The frequency values for 
fibre -orientation angles 15 ° and 30° lie in between 1he solid line and the 
dotted line for a particular modulus ratio, but are not shown. 

The variation of andm^^ with mo.dulus ratio for c/a = 0.0 

is shown in Figure 2.22 for a square, simply-supported plate. 'Fot all 
fibre -orientations and increase with modulus ratio* However, 
and cD^-j behave in exactly opposite mamer; for ary modulus ratio, ! 

increases with increase in fibre -orientation angle whereas decreases, i 

I ! 

This trend is found to be true for all c/a values, 
d) The effect of cut-out on the mode shape-. 

The oigcnvcctors, normalized with respect to the irass matrix have 
been obtained in all the cases. It is found that symmetric modes_ are mor^ 
influenced due to the presence of a cut-out. The coupling between modes , 
for a symmetric mode is found to be more pronounced at higher x cut-out i 
ratios. This seems logical as the cut-outs considered are symmetrically, located. 
The antisymmetric modes and a corbination of symmetric and antisymmetric 
modes are not influenced much. Table 2.6 gives a comparison of frequency 
values (for a symmetric mode) and the normalized eigenvectors for all 
modulus ratios in the case of 6 = 0°. The influence of cut-out can be seen 
on the eigenvectors. 
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It has also been pointed out by Ivlaiiabaliraja et*al. CssJ'^hat the 
cut-out has an effect of coupling the different waveforms in an 
otherv^ise solid shell. Ihis has been observed to be true in the present 
case of plates vd.th cut-outs also and therefore, the identification of 
eigenvalues from the largest amplitude coefficient in the eigenvector 
becomes difficult. However, it is possible to identify the eigenvalues 
from a comparison of corresponding eigenvectors and finding the relative 
values of anplitude co -efficients in these eigenvectors. 

2 , 4. 3 Clamped -Clamped Plates 

a) Effect of cut-out on the natural frequency 

The frequencies and corresponding eigenvectors have been obtained in 
all the cases considered, A verification of results with available results 
is made for an isotropic solid plate* Table 2*7 gives a comparison of the 
frequencies and corresponding eigenvectors for a clamped-clamped isotropic 
solid plate. It is seen that the results obtained in the present investi- 
gation compare veiy well with Young^s £832 results , 

Table 2.8 gives a comparison of the non-dimensionalized froquenci® 
of an isotropic plate for cut-out parameters equal to 0.0 and 0 . 5 . It is 
observed that the value of the &ndamental frequency for a plate without 
cut-out obtained in the present study is the same as obtained in 18 43 * 
However, the value obtained in C52l is even less than the lower bound 
reported in 11841. This indicates the improved accuracy of the present I 
method over other methods. It is also seen that the values of frequencies | 
for c/a = 0*5 reported in P" lower than the values obtained xn the 



Table 2.7 

Comparison of frequencies and modes of vibration of an isotropic clamped-clamped 

square plate; c/a = 0.0. 
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Table 2.8 

Comparison of non-dimensionalized frequency for a clamped-clamped, 
square isotropic plate; c/a = 0.0 and c/a = 0.5 


Cut-out parameter 

Source 

0.0 0.5 



Mode 1 

Mode 1 

Mode 2 

Mode 3 

Mode 4 Mode 5 

Kristiansen 


60.80 




CsO 

Paramsivain 

34.85 

62.40 

68.55 

98.52 

143.95 

C 52 II 

Agaiwal 

35.7 

69.10 

77.20 

105.0 

— — 

Casq 

leissa 

C843 

35.985 

(lower bound - 
35.976) 

■ 




Present 

35.985 

65.581 

79.619 

108.91 

117.95 190.83 
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present study; this may be due to the fact that a coarse mesh size has 
been taken by Paramsivam. A comparison is also made with the availa,ble 
experimaital results j s^Jand good agreement is £)und. 

The. variation of non-dimensional frequencies with c/a for a clamped- 
clamped isotropic plate is sho\^ in Figure 2.23. It is interesting to note 
that the third and fourth modes as well as ± the fifth and sixth modes 
intercharge. This has also been observed in []86 ^ foi* a- rectangular iso- 
tropic plate with a central square cut-out. Results obtained in Da are 
lower than those obtained in the present study. The values given by 
Kristiansen and So e del compare well for smaller values of c/a but are 
lov^ for hi^er values of c/a. This may be because of the reason that they 
have used Hayl eights method. The experimental values given in Co 53 are 
also compared for the first three natural frequencies. 

The fundamental frequency, unlike in the case of a simply -supported 
isotropic plate, does not decrease for small cut-outs but increases conti- 
nuously as the cut-out size is increased. The behaviour of frequencies 
^12 ^ " ^21^’ ^22 ^^^*^^23 ^ ~ ^32^ found to be similar. They first 

decrease with the increase in cut-out size but again increase as the 
cut-out size is increased. This has also been observed for a simply- 
supported plate. It was also seen that the fourth mode and the fifth mode 
are identical for a simply -supported solid isotropic plate. This is not j 

found to be true for a clamped- clamped plate. The reason is that a closed j 

I 

form solution for a clamped- clamped isotropic plate is not known and the : 
approximation can be assumed to gi-ire rise to the fourth mode even when 
there is no cut-out in the plate. 
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In Figures 2.24 and 2.25} the variation of with c/a is shovsn for 
balanced bidirectional composite plates in the case of 6= 0° and 6 = 45° 
respectively. It is fcund that and increase with c/a for = 0° 

and e = 450 as well. For 6 = 0° and 8= 45°, '^^ 2 ’ ”^22’ '^13 ^23 

decrease with incareasing values of c/a and then increase again for higher 
values of c/a. It is i^eresting to note that modes corresponding to 
^22 ^**13 those corresponding to and interchange for aH 

fibre-orientation angles. The behaviour of balanced bidirectional composite 
plates is veiy similar to that of anisotropic plate. 

The variation of m with c/a for unidirectional glass-epoxy plates 
is shown in Figures 2.26 and 2.27 for6= 0° and 6 = 45° respectively. Here 
also, like balanced bidirectional plates, and increase with 
increasing values of c/a. For 8 = 0°, modes corresponding to and 
interchange at a smaller value of c/a. As the fib re -orient at ion angle is 
increased they separate from each other at smaller values of c/a and 
interchange at a higher value of c/a (Figure 2.2?). The modes correspond- 
ing to cogj “j-] not interchange when 6=0° (Figure 2.26) but as the 
fibre-orientation argle is increased they cane closer and separate from 
each other after intercharging once with the increase in c/a ratios 
(Figure 2.27). As the fibre -orient at ion an^e is increased, the frequency j 
curves for ^nd and those for ^nd come closer and coincide j 

for 0 = 45° . 

The variation of frequencies with cut-out ratio is shown for a 
uni directionally reinforced boron-epoxy plate in Figures 2.28 and 2.29 for 
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9=0° and 0 = 45° j resp cjctively. The frequencies sJ3.d 

increase with c/a ratios for both 6=0° and 0 = 45°. The behaviour of 
frequency is, however, different; it increases with c/a ratios when 

0=0° (Figire 2.28) whereas it is feund to decrease with inci-easing 
values of c/a and then to increase at higher values of c/a when 6 = 45° 

(Figure 2.29). The behaviour of frequencies “^2’ ^21’ ^22^ ^14 ^23 

are found similar both fcr 6 = 0° and 6= 45°. Ihe mode corresponding 

to is found to intercharge with those corresponding to ^ 2 -\^ ^22 

v\?hen0= 0°. However, ihere is no interchange of modes except for 

those corresponding to and in the case of 6 = 45°. The frequency 

curves for and oj., and those for lOo-r and a)__ come closer with the 
i2 21 23 32 

increase in fibre -orientation angle and they coincide for 0 = 45°. 


Fo r unidirectional graphite -epoxy plates, 


the variation of to with c/a 


for 6 = 0° is sbown in Figure 2.30. The frequencies 

(correspondiag to symmetric modes) increase with c/a. The behaviour of 


frequencies and “^22 is similar to that found for balanced- 

bidirectional, unidirectional ^ass-epoxy and unidirectional boion-epoxy 


plates. The mode corresponding to ^22 interchanges with that corres- 
ponding to and then with modes corresponding to <^21 > ^^14 '*'13“ 

0 = 45°, the variation of m with c/a for imidirectional graphite-epoxy 


plates is shown in Figure 2.31. Here also, as for 6= 0°, increases 
with c/a. The frequencies ^^2? ^22 found to decrease with 

increase in c/a values; the decrease is significant in compared to 
that in aaccase of 6 = 0°. These frequencies, however, increase with hi^er 
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values of c/a. It is interesting to note that there is no_ interchange 
of first four modes when. 8 = 45°. It is observed that for clamped-clamped 
plates ?/ith contra! square cut-outs, the frequencies for each mode except 
the first vary with cut-out size in different ways and do not uniformly 
increase or decrease, however, as in the case of simply -supported plates 
some similarity is observed in their behaviour, Ihe fundamental fre- 
quency in all the cases considered is found to increase with increasing 
c/a values. This suggests that cut-outs make plates stiffen for clauped- 
clamped edge conditions. Frequencies “^ 2 ’ ^ 21 ^ '*^22 ^23 found to 

first decrease with increase in cut-out size for all fibre-orientation 
angles and modulus ratios . They, however, increase for larger cut-out 
ratios. The frequency curves for ^21 those for and 

are found to come closer as the fibre-orientation angle is increased and 
these coincide when © = 4 . 5 ° for all the modulus ratios (this behaviour is 
found in an isotropic plate). The modes corresponding to and 

are not identical in this case as was seen in a simply -supported solid 
plate for fibre -orientation angle equal to 45°. 

b) The Effect of Pib re -Orientation Angle on the Natural Frequency 

The influence of fibre-orientation angle has been shown only 
for for different values of c/a ratios. The variation of 

with 0 for different c/a values is showi in Figures 2.32 to 2,35 
for various composite plates. In Figure 2.32, which shows the variation of 
with 0 for balanced bidirectional plates, it is observed that the rate 
of decrease of frequency increases with cut-out parameters. This suggests 
that the larger cut-outs in this case tend to make the plate less stiff 
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when 0 is increased. In the case of unidiroctional glrsss-epozy plates 

(Figure 2,33)5 <iecreases with increasing 0 for medium cut-out sizes. 

As c/a increases 5 increases continuously v/ith increasing values of 

6 (c/a = 0,65 c/a = O.7). The behaviour of with 0 for a unidirectional 

boron-epoxy plate is found similar to a unidirectional gla.ss-epoxy plate 

and is shown in Figure 2 , 34 . In Figure 2,3 5^ which is a plot of vs, 

® for unidirectional graphite-epoxy plates j, it is seen that for small 

values of c/a (upto 0 . 2 )^ decreases with increasing 0 , For values of 

c/a lying between 0,3 and 0,45 it first increases upto a value of 0 equal 

to 22 , 5 ^ and then decreases. As c/a is further increased, . is found 

^ ^11 

to increase with 0 , 

It is seen that the variation of the fundamental frequency with 
0 for small values of c/a (upto 0,.2) is similar for all composite plates; 
it is foimd to decrease with incroa.se in the fibre -orientation an^e. Thus 
increase in the fibre-orientation an^e makes clamped— oLampod plates with 
small cut-outs less stiff. At higher c/a values ^ the fundamental frequency 
except in the case of balanced bidirectional plates increases with the 
fibre-orientation an^o, Maurizi and Laura DsD also observed the 
decrease in the fundamental frequency with fibre -orientation angle for a 
square 5 clanped-clamped boron-epoxy plate with no cut-out, 

c) The effect of modulus ratio on the natural frequency 

The influence of modulus ratio on the fundamental frequency of a 
plate with a cut-out is shown for two fibre-orientation angles - 0^ and 
45 *^. For different modulus ratios ^ the variation of with c/a is shown in 
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Figure 2,36 in the case of 0= 0° and 6= 45° • It is found that increases 
with c/a for all aodalus ratios. The frequency curves for 9=0° and 
6 = 450 for a given modulus ratio, except in the case of balanced- 
bidirectional plates, are found to intersect with each other. This indicates 
the increase in stiffness, at higher c/a values as the fibre-orientation 
angle is increased, for plates with ratios of principal Young's modulus 
equal to 3, 10 and 40. It is also observed that the frequency values for 
an isotropic plate are higher as compared to those for a baladiced-bidirectional 
plato in the case of 6 = 0° and 9 = 45°. 

The variation of 03^^ andu^g wi'th modulus , ratio for c/a = 0,0 is 
shown in Figuare 2,37. For all fibre -orientetion angles, and 
increase with modulus ratio. But for a given modulus ratio, and 

* behave in exactly opposite manner.: is found to decrease with increase ^ 

in fibre-orientation angle whereas is found to increase v^ith incareaso 
in fibre-orientation angle. 

d) The effect of cut-out on the mode shape 

The eigenvectors noimalizod with respect to the mass imtrix have been 
obtained in all the cases considered. It is found that symmetric modes are 
more influenced due to the presence of a cut-out. This is because of the 
symmetric location of the cut-out. Table 2.9 gives a comparison of 
frequency values for a symmetric mode and the noimalizcd eigenvectors for 
all modulus ratios in case of 0 = 0°, 
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CHAPT3E 3 


VIBHATIONS OF' GOtlPOSIfS PLATES Y/ITH HOLES 


3.1 Introduction 

A detailed study of the vibrations of plates with cut-outs was done 
in Chapter 2. It was noted that the theoretical results for the free 
vibrations of pla.tes mth cut-outs could be compared only with the 
available theoretical results for isotropic plates and a few experimental 
results available on clamped-clamped isotropic plates; practically there 
are no results availa.ble on vibrations of composite plates with cut-outs. 
The theoretical study on the vibrations of composite plates with holes 
has been done with a view to verify the computed results with experimental 
results. It is also easier to make holes in plates than to make cut-outs 
as the latter involve practical difficulties, especially in composite 
plates. In this chapter, the theoretical results are given for sq.uare 5 
simply-supported and clamped -clamped composite plates with centrally 
located circular holes. Ho results are given here for balanced bidirec- 
tional ^ass-epoxjr plates as it was fotmd (in Chapter 2) that the behaviour 
of balanced bidirectional ^ass-epoxy plates is very much similar to that 
of an isotropic plate. A comparison of results has also been made for 
isotropic plates with available results and the experimental results 
obtained in the present study; the experimental study is discussed in 


detail in Chapter 4 
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3»2 Method of Solution 


It was seen in the vibra.tions of plates T>;ith cut-outs that the 
integrals necessary for the evaluation of the coefficients 

^mn’ ^im *^3110 found in dosed form. There is no change 

in the expressions for and and they are given by equations (2.19, 

2.21) and (2.30, 2.32) for a simply -supported and a clamped-claEped plate 

respectively. Hov/ever, the integrals necessary in the computation of 
i j 1 j 

obtained in closed form and hence a numerical 
integration scheme is necessaiy* Of the various schemes ^ the sixteen 

quadrature formula has been found to yield sufficiently accu- 
rate results and has been used in the evaluation of 6 ^^ and The 

mn mn 

coefficients and arc given below for simply -supported and 

clamp ed— clamped boundary conditions. 

3 . 2.1 Simply-Supported Plates 

Ihe expressions for 3^ and for a simply-supported plate are.^ 



== 

phab 

%1 



(3.1) 
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where , 

jj sin (inirX) sin (irrX) sin (n ttY) sin (3 ttY) (Ev dY 

jj cos (m irX) cos (iirx) cos (niTY) cos (3 ttY) dil dY 

// sin (mTi’x) cos (i irx) sin (n ifY) cos (jirY) dS dY 

// cos (mTfX) sin (i ttx) cos (n tty) sin (j''^Y) dX dY (3.3) 

Where Att is the area of the hole . 
xl 

X = x/a, Y = y/b . 

D^^'s are the elements of the plate bending stiffness matrix and 
are related to the orthotropic elastic constants of the layer and the 
fibre orientation angle 0 , of the layer with respect to the geometrical 
axis of the plate. 

3.2.2 Clamped-Clamped Plates 

The expressions for 3^^ and for a clamped -clamped plate are 

mn mn 

6^^ = P hab // X X. Y Y. dX dY 

mn ' mini 


f? 


"HI 


^H2 


■H3 


P, 


H4 


(3.4) 
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where 5 


a 


,13 

mn. 


D22^b 


a 1- 22 

2 ®66 


^ ^ * V * *’ ^4H 


22 


D 


16 


+ 4E ■=: ^CTT + 2R ^ i- C.-CJ ^VTi^ 


^22 5H 


D22 "eH 'IE' 


, 2 R^ (ji 4 . 1? ) 

^ ^ ©22 ^ 8H "9H^ 


(3.5) 


■1H 


- rr X X Y Y. <ax dY 
~ JJ m.XX ijXX n 3 


A 


■p =/’fxX Y Y.dXdY 
^2H m i,XX n,YY 3 


P -f/x XYY dXdY 
^3H “ m^XX i n "^jjYY 

Ajj 


F 


4H 




dX dY 




S' = ff X ^ X. ^ Y „Y.„dXdY 
5H *'4 ia>X; i,X n,Y 3 jY 


^6H " ^m,XX ^i,X ^3,Y 


dX dY 




p - (f y: X Y YdXdY 

^7H “ m,X i,XX n,Y j 

^8H " \,X ^n,YY ^3,X 

H 


% = IJ \,X "i ^n,T ^j.YT “ ® 

H 


(3.6) 



3.2.3 Evaluation of and 

nn Em 

The integrals given in equations (3.3) and (3.6) for the evaluation 
sinply-supported and clanpe d-clampe d boundary- 
conditions in the case of a cut-out have been expressed in closed fom. 
However 5 for a hole there is no such sinplification and a numerical 
integration scheme is necessaiy to evaluate these integrals. A sixteen- 
point Gauss quadrature formula has been used for this purpose. The 
sixteen-point formula has been compared v/ith t-wo nine-point Gauss quadrature 
formulae to evaluate a standard integral for which the exact solution is 
known and it was found to give more accurate results than the nine-point 
formulae. Table 3.1 gives the three numerical integration schemes with 
numerical quadrature points and corresponding weights. Table 3.2 gives 
the values and the coniputation time required on an IBM 7044 computer with 
these numerical in-tegration schemes to evaluate the standard integral. 

The percentage error using the sixteen-point formula in the evaluation of 
the standard integral is found to bo less than 1.2 percent and this error 
is considerably less than that obtained by using the nine-point formulae. 
Therefore, the integrals in equations (3.3) and (3»6) have been evaluated 
using a sixteen-point formula even though the total computation time for 
the evaluation of the standard integral is more than that required when 
nine-point formulae are used. The integrals evaluated once for a parti- 
cular hole size, are stored in the memory and used whenever needed. The 
use of the numerical Integration scheme increases the computation time 
required for the determination of natural frequencies and corresponding 
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Table 3.1 


numerical quadrature points and corresponding weights for different 

Gauss -quadrature formula 

n 

// f(x,y) di. = A I W. f(x., y^); n = 9 or 16 . 

A i=1 


Gauss Quadrature 
Ponrula 


line -Point (l) 


nine -Point (2) 


Sixteen-Point 


i 


^i 

‘‘i 

1 

0.0, 0.0 


0.25 

2 

0.92387953, 

0.38268343 

0.09375 

3 

0.92387953, 

-0.38268343 

0.09375 

4 

-0.92387953, 

0.38268343 

0.09375 

5 

-0.92387953, 

-0.38268343 

0.09375 

6 

0.38268343, 

0.92387953 

0.09375 

7 

0.38268343, 

-0.92387953 

0.09375 

8 

-0.38268343, 

0.92387953 

0.09375 

9 

-0.38268343, 

-0.92387953 

0.09375 

1 

0.0, 0.0 


0.16666667 

2 

0.70710678, 

0.0 

0.16666667 

3 

-0.70710678, 

0.0 

0.16666667 

4 

0.0, 0.70710673 

0.16666667 

5 

0.0, -0.70710678 

0.16666667 

6 

0.70710678, 

0.70710678 

0.04166666 

7 

0.70710678, 

-0.70710678 

0.04166666 

8 

-0.70710678, 

-0.70710678 

0. 04166666 

9 

-0.70710678, 

■ 0.70710678 

0.04166666 

1 

0.42470820, 

0.17591990 

0.0625 

2 

0.17591990, 

0.42470S20 

0.0625 

3 

-0.17591990, 

0.42470820 

0.0625 

4 

-0.42470820, 

0.17591990 

0.0625 

5 

-0.42470820, 

-0.17591990 

0.0625 

6 

-0.17591990, 

-0.42470820 

0.0625 

7 

0.17591990, 

-0.42470820 

0.0625 

8 

0.42470820, 

-0.17591990 

0.0625 

9 

0.82047322, 

0.33985115 

0.0625 

10 

0.33985115, 

0.82047322 

0.0625 

11 

-0.33985115, 

0.82047322 

0.0625 

12 

-0.82047322, 

0.33985115 

0.0625 

13 

-0.82047322, 

-0.33985115 

0.0625 

14 

-0.33985115, 

-0.82047322 

0.0625 

15 

0.33985115, 

-0.82047322 

0.0625 

16 

0.82047322, 

-0.33985115 

0.0625 
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Table 3.2 

Comparison of results and computation time to evaluate a standard 
integral using various numerical integration schemes. 



I = 

11/, 

/ - 

0 0 

dx dy 


Integration 

Exact value 

Computed 

Percentage 

Computation time 

S cheme 

of I 

Value 0 f I 

error 

on an IBL' 7044 

Computers 

in milliseconite 

1 

2.0943951 

2.3561944 

12.6 

616 

2 

2,0943951 

2.0045597 

- 4.3 

616 


3 


2.0943951 


2.1170794 


1.13 


716 
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Gigenvcctors. For example^ the total computation time required for the 
deteriDination of four natural frequencies and corrospondiiig eigenvectors 
considering only symoatric inodes is approxiiBately nine minutes for a cut- 
out and Gighte‘en ninutes for a hole and all fibre-orientation angles 
(o® 5 15^? 30° j 45®) in the case of clamped-claiiped plates. 

3*4 Results and Discussion 

Natural frequencies and corresponding eigenvectors have been 
computed for square^ simply-supported and clanped-clamped plates with 
centrally located circular holes for different fibre-orientation angles 
in steps of 15®# Results for balanced bidirectional plates have not been 
given as the behaviour of natural frequencies for those plates is 
similar to that found in the case of an isotrapic plate. Ratios of 
principal Young ^ s moduli o'f 3^10 and 40 have been used in the computation. 
Results are given in terms of non. dimens icnal frequency OJ for fibre- 
orientation angles 0® and 45® only* A comparison of frequencies and corres 
ponding eigenvectors for a particular mode is made between cut-outs and 
holes to study the influence of other modes, Ihe oDn'\^rgence of the 
results has boon studied and a relative comp axis cn of the frequency 
values led to the decision to compute natural frequencies taking six'been 
frequency equations and different combination of modes for ail the cases 

considered. It should be noted that twenty-five frequency equations y/ere ! 

) 

taken in the computation of natural frequencies for plates with cut-outs,; 
The number of equations has been reduced for plates with holes to make a 
coii$>roiiiiso between the computation time and the accuracy involved. Also^ 
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computational difficulties like non-convergence of frequency values ^ 
increa.se in round-off errors were encomtered at larger hole sizes and 
therefore results could be computed for the hole sise parameter vai^^ing 
from 0.0 to 0.5 only except in the case of fandariental frequency# 

3.4 #1 Simply -Supported Plates 

a) Effect of Hole on the Natural Frequency 

Table 3*3 gives a comparison of the non-dimensional frequencies 
in the case of a square^ isotropic platewith a centrally loca.ted circular 
hole v/ith available results. The lesults given by Kumax O-bH 
terms of frequency ratios have been converted in the non-dimensional 
form for cori:q)aiiso n purposes. There is a good agreement between 
Kumai^s results and the results cemputed in the present case. However^ 
frequency values obtained by Kumai are lov^er than the values obtained 
in the present study. Kumai computed frequencies by the superposition 
of the solutions of circular ring plate satisfying the bounda.iy conditions 
at several points along the edges of the square plate. This cpp 2 X)xiiaation 
allows some residual deflections at some portion of the outer boundaiy of 
the plate vthich can be assumed to make the plate less stiff and hence 
low frequency values can be expected. Joga Rao and Picket [3-83 used 
Rayleigh-Ritz method and modified the admissible functions to take into " 
account the appropriate singularities within the holes. A comparison | 
with the present theoretical and experimental results shov/ the inaccu- 
racy of the result obtained in Reference [^483 sven after making 
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corrections for the singiilaiitios in the holes and thorcforcy it is 
doubtful whether the method suggested there can give accurate results. 

The variation of frequency with hole parameter ^ d/a for an 
isotropic 'Plate is sliewn in Figure 3 . 1 . The frequencies have boon 
subscripted (11^ 123... etc.) on the basis of the largest amplitude 
coefficients in the corresponding normal i 2ed eigenvectors, A good 
agreement is obtained v/ith the experimental results. The presv;nce of 
hole in a plate gives rise to a new mode similar to the case of a plate 
with a cut-out. For the d/a ratios considered 9 there is no interchar^e 
of modes. 

In Figures 3.2 and 3.3 the variation of o) with d/a is shown for 
unidirectional glass-cpoxy plates in the case of 6=0^ and 6 = 45 ^ res- 
pectively, For 6 = decreases veiy little with the incox^ase in 

d/a values and then increases as the hole size is increased , There is no 
decrease in with d/a for 6 = 450, it is scon that as the flbrc- 

oilentation angle is increased, frequency curves for ^^2 ^21 

closer and they coincide v/hen 9 = 45°? this is also true for frequency 
curves and The rrodcs corresponding to and seem to 

have a tendency to intercharge at d/a = 0 , 5 . 

The variation of ^ with d/a is shov/n in Figures 3.4 and 3*5 for 
unidirectional boron-epoxy plates. It is seen that for 9 = 0 ^ (Figure 3 * 4 )^ 
all frequencies except decrease with increase in d/a values and at 

higher d/a values only and ^2^ increase. The modes corresponding to 
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and 05^^ interchan;® and those corresponding to and ^22 ^ 

tendency to interchange. For 6 = 45°? there is no docroase in 
(Figure 3.5). h^)dcs corresponding to Wgg found to inter-- 

chan^ at higher d/a ratios. 

For unidircctionally reinforced graphi to -epoxy plates, Figures 3.6 
and 3.7 shov/ the variation of w with d/a for 6 = 0 ° and 6= 45°, vYhen 
0=0°, and decrease with increasing d/a ratios and the rate 

of decrease is more as compared to isotropic, unidirectional glass-epoxy 
and unidirectional boion-epoxy plates. However, for 6 = 45°, ^ 

increases vhLth increaso in the hole size. The frequenev curves for 
and come closer as the fibre -orientation angle is increased and these 
coincide when 6 = 45°. This trend was also found in unidirectional 
glass-epoxy and unidirectional boron-epoxy plates. 

b) Effect of Fibre -Orientation Angle on the natural Frequency 

The effect of f ibro-orientation angle on the frequencies and 

“^2 been shown in Figures 3.8 - 3.10 for selected values of d/a ratios. 
It is scon that in all the cases and increase with the fibre- 

orientation angle for all d/a ratios. As the hole size is increased, the 
rate of increase of increases. For Tonidircctional glass-epoxy plates 

(Figure 3.8), corresponding to d/a = 0.3 is lower than that correspond 

ing to d/a = 0.0 at 6 = 0° and the curves for corresponding to these 

two values of d/a intersect as 0 is increased. In Figures 5.9 and 3.10, 
a similar behaviour is observed for unidirectional boron-epoxy and 
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graphite-epoxs^ plates, The frequency found to increase with 0 

for all hole parameters similar to also observed from 

Figures 3.8 to 3.10 that fora given value of 0, the frequency 
decreases with increase in d/a values whereas this behaviour was observed 
in case of .only for smll values of 0 and valuesof d/a upto 0.3. 

c) The Effect of Ibdulus Ratio on the Ratuial Frequency 

The influence of modulus ratio on the fundanental frequency of a 
plate with a hole is shown for two fibre-orientation an^es - 0° and 45°. 
For different modulus ratios, the variation of with d/a is shown in 

Figure 3.11 in the case of 0 = 0° and 6 = 45°. It is observed that for 
isotropic plates, decreases for small d/a ratios vdicreas for an 

orthotropic plate (6 =0°), the decrease in occurs upto larger 

values of d/a. The rate of decrease increases with increasing modulus 
ratios. The variation of with d/a for different modulus ratios is 

shown by a dotted line when 0 = 45°. The fundamental frequency is foxmd ; 
to increase for all modulus ratios with d/a. 

d) The Effect of Holes on the Mode Shape 

For plates with holes also, the eigenvectors normalized with 
respect to the mass matrix have been obtained. A comparison of frequency 
values and corresponding eigonvoctors has been made (Table 5.4) for the 
same hole and cut-out parameters for a symmetric mode. The cc 2 i"cribution 
of other modes is found significant for a cut-out as well as for a hole 


in the plate. 
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6, degrees 


Fig. 3. 8 w vs. e for different hole parameters for 
a simply-supported plate;R=10. 



0, degrees 

Fig.3.9 55 vs.© for different hole parameters for 

a simply-supported plate;R=1*0. 



1 /iD^2 



d/a^' 

Fig. 3.11 05 ^^ vs. d/a for different modulus ratios 

for a simply-supported ptate;R=1'0. 



116 


rr\ cvj cr\ ltx trs 
lA tA 'T- CA VD CO '==-;r VO 

VD i>- CO O CA I>- CM <“* 

VD 'vh T- Cvi O T- T~ O 


o o CO CO 

CO OA <- CkI IA- CO 

^ A- VO CM O CJ eg 

o o o o o o o 


V- o o o o o o 

I 1 i I 


o o o o 

I I 


O K\ 'T-' CA [>. VO o 

O VO T- v- CA LA eg 

['A O tA O OA CA CA 

T~" LA CM CO CO "T" T~“ 


VO tA t-A f>- LA O LA 

CA CM LA CM t — — VD 

T- V- rvA LA - LA O 

o o o o o o o 


G O O O O O 

! I ! 


o o o o o o o 

I { I 


32 J 

} i 

I 

A- 

' CM 

LA 

V- 

o 

I A 

G 

A- 

T” 

^.A 


LA 

LA 

VD 

CA 

VQ 

LA 


! I 

11 1 

O 

■ 

VO 

GA 

o 

\ — 

LA 

LA 

!>“ 

CM 

o 

LA 

O 

A^ 

CO 

JA- 

CM 


1 1 

! 

• 

' LA 

tA 

eg 


LA 

oo 

LA 


O 

o 

CD 



o 

CM 

LA 

<N 

I { 

Cl 

o 

: CA 

G 

LA 

CM 

■'i — 

o 

CM 


T — 

o 

O 

o 

CO 

* 

o 

O 

O 

# 

o 1 

I ! 

\l 

o 


m 

• 

« 

« 

• 

# 

« 

# 

• 

• 

♦ 


• 

O 1 

lOi 

1 

V- 

* O 

T— 

o 

o 

o 

o 

o 

O 

o 

o 

o 

o 

o 

{ 

CO 

o 

O 

1 

fv— I 

1 


^ I 


1 


i 



I 


1 

I 



1 

11 : 

1 1 

lA! 


' 















a> J 

CM! 1 

• 1 

O! 


i 

















!C£j 1 1 

I 

CM 

^ V- 

VO 

T— 

GA 

LA 

(A- 

lA- 

o 

LA 

tA 

o 

eg 

\D 


CM 

(A- 

1 

2 1 

I! 1 

T~ 

■ 

LA 

lA- 

CA 

exo 

CA 

CA 

OA 

LA 

O 

' g** 

LA 

LA 


LA- 

T— 


T— i ! 

1 

05! 

« 

! i>- 

'-■g- 

NA 

O 

CA 

T— 

£A- 

CO 

O 

VO 

OJ 

VO 

eg 

CO 

VO 



i-M * * 

CO 


O 

— 

CO 

o 

A 

lA- 

1A~ 


eg 

o 

CM 

lA 

CM 

• 

T — 

o 

0 ffv ^ 

! I 


LA 

1 ♦ 

• 

• 

• 

» 

# 

« 

* 


« 

• 

• 

* 

• 


1 I 

TO 1 


i T“ 

^ — 

V— 

o 

o 

O 

O 

O 

o 

o 

o 

CO 

I 

O 

O 

o 

o 


! I 

1 


! 5 



I 



1 


1 




1 

! 

? 

1 1 

f 1 

lAl 


4 















* 4 

ii I 

i i 

1 i 

♦ t 

OI 


* 

















1 [ 

1 

lA- 


VO 

LA 

eg 

{A- 

O 

VO 

CA 

LA 


"h 

o 

CO 

CD 

00 

OA 

o ’ 

I 1 

il { 

CA 

1 ^ 

CA 

LA 

T“ 

VD 

LA 

VO 

'•'A 


6 

VO 

CM 

LA 

OA 

-■A 

CA 

1 

i i 

I f 

1 

o5 1 

•s — 


OA 


CD 

VO 


o 

CM 

OA 

LA 

LA 

VD 

CO 

NA 

O 

— 

05 

1 I 

m 

] VO 

LA 

■«“ *- 
"4 

r— 

<- 

o 

CM 

— 

v- 

o 

O 

O 

O 

O 

O 

o 

# 

1 

1 I 

\l 

CO 


• 

« 

« 

• 

« 

# 

m 

• 

• 

• 

• 

• 

* 

• 


I { 

O 1 

CA 

^ O 

V“ 

o 

o 

o 

o 

O 

O 

o 

o 

o 

o 

o 

o 

o 

o 

-a 

ILAI 

I 


! * 


1 


I 

1 


r 


I 

I 


i 



I 


S 1 

ITst 


1 

i 
















I I 

I r oi 



0 

1 ! 

1 

LA 

CA 

VO 

o 

CA 

lA- 

O 

CM 

o 

VD 

LA 

OA 

VD 

A- 

LA 


O 

05 

i — 1 


II i 

CM 

G 

lA- 

CA 

CM 




LA 

LA 

O 

•4~ 

— 

VD 

CA 

CA 

■\ — 


o 


f 

• 

CM 

T- 

T— 

[A- 

O 

A- 

LA 


CM 

o 

o 

CM 


LA 

T— 

o 

m 



Co * 


LA 

T— 

o 

O 


O 

O 

O 

o 

o 

O 

<D 

O 

o 

o 

CD 




o 

• 

m 

• 

• 

• 

« 

W 

» 

0^ 

• 

* 

• 

• 

• 

# 

•> 


o5 

J 1 

«o 1 

T- 

O 


o 

O 

o 

O 

o 

O 

o 

o 

o 

O 

O 

o 

o 

o 

P_j 


! ! 

s 


I 


i 

1 

I 


1 





1 


I 


I 

05 


1 { 

LAI 

















1> 

_P 

1 ! 

« { 


















l>> 

•H 

I ! 

Ol 



















1 i 

I 

-A 


CA 

LA 

LA 

CA 

T- 

CA 

CO 

LA 

LA 


o 

LA 

CO 

o 

o 



1 1 

il I 

CM 

lA- 

LA 

VO 

T— 

LA 

CM 

LA 

— 

VD 

LA 

LA 

LA 

T— 


CA 

CA 




1 

* 

T- 

-g- 

CM 

CA 


O 

'K — 

— 

CM 


VO 

tA 

m 

— 


CA 

0 


* OI 

C^l 

VO 

LA 

o 

CM 

O 

o 

LA 

o 

o 

CM 

O 

O 

o 

o 

o 

o 

O 



I-HI 



• 

« 

* 

• 

» 

• 

• 

« 

• 


m 


• 

« 

* 

« 

QTl 


IP-i 

O i 


V- 

T — 

O 

O 

T— 

o 

o 

o 

O 

o 

a 

o 

o 

o 

o 

o 

0 


I Ol 

1 


1 


I 





1 

1 



1 


I 

i 


u 

i 

11^1 

LAI 

• { 




















I OI 

O J 

4 



















1 mi 

t 

CA ' 

















O 

i 

iMi 

11 1 

£>- J 


















1 


1 

* 1 



















It- I 

Ci;5 1 


CA 

o 

o 

CA 

o 

e- 



9 


o 

VO 

CA 


VO 


O 


I t 

"■^1 

''■i j 

■s; — 

T— 

Tu- 

VD 

T— 

VD 

o 

LA~ 

'Or 

o 


o 

VO 

A- 

o 

00 

m 


t 1 

rai 


CA 

LA 

"V — 

O 

LA 

lA- 

T— 

o 

T— 

T- 

O 

t: — 

o 

o 

•t — 

o 

•H 


j 



OA 

CO 

o 

O 

CO 

LA 

o 

o 

O 

o 

o 

o 

o 

o 

o 

o 

U 

5 



LA 1 

• 

« 

« 

* 

• 

« 

• 

9 

* 

* 

• 

♦ 

• 


*' 

• 

o5 


i 



O 

o 

o 

O 

o 

O 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

S 

o 

1 

I ' 


1^ 

I 






I 

! 


I 

I 



! 


1 

o 

j 

f 


•r-j 

T- 

LA 

LA 

tA- 


LA 

LA 

lA- 

T— 

LA 

LA 

^ LA- 


lA 

lA 

^ A 



i 


'H* 


■M — 

T— 

T— 

LA 

LA 

LA 

CA 

LA 

LA 

LA 

LA 

, £A- 

A^ 

A- 

A 





<r* ' 

! <4 



<4 


- J 


<4 

<E[ 


<3^ 

cej 

<s:l 


<4 




117 


3.4 <.2 Clampcd-Clainped Plates 

a) Effect of Hole on the natural Preg^uenej 

A conpa.rison of the non-dimens ional frequencies for a few node has 
been made in fable 3.5* The results obtained by Anderson i3o3 ^7 using 
finite element method compare well with the present results. 

The variation of o) with d/a for a clasped-clanped isotropic 
plate is Bhowi in Figure 3.12. The third and fourth modes are not found 
to interchange unlike in the case of clanped-clamped isotropic plates 
with cut-outs (Figure 2 . 22 ) whereas the 5 th and 6 th modes are found to 
interchange with each other* Eesults obtained by Anderson {3oII compare 
well with the fundamental f requency whereas those obteined by Kumai [ 4 6 ^] 
are lower. It was mentioned in Section 3.4.1 that the approximation of 
super -posit ion of the solutions of circular ring plate allov\^s some 
residual deflections at some portion of the outer boundary of plate which 
in facty makes the plate less stiff. Therefore ^ lower frequency ■values 
are obtained by Kumai. The experimental results could not be obtained 
for the third aid the fourth modes? however^ the experimental results for 
other modes compare well with the theoretical results. 

In Figures 3.13 and 3.14, the vaiiation of w with d/a is shewn 

for unidirectional glass-epoecy composite plates in the case of 8 = 0 ® and 

e = 450 respectively. The frequencies w.. and are found to 

n 31 

increase with d/a ratios. For 9 = 0^, modes corresponding to and ^ 

interebarge at a sebII value of d/a. As the fibre-oiuentation an^o is 
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Table 3-5 


Comparison of non-dimensional flrequency for a square clamped-danped 
isotropic plate with a central ■ circular hole. 


d/a 


Present 

Anderson 

■first Ko de 

Fourth Mode 

First Mode 

Fourth Mode 

0.0 

55.985, 

131.589 

35.375 

130.540 

0.125 

36.321 

129.881 

35.138 

128.637 

0.25 

37.732 

128.200 

37.600 

125.464 

0.375 

43.752 

125.600 

44.240 

122.020 

0.5 

55.610 

125.212 

57 .844 

124.308 
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increased they separate fron each other* Modes corresponding to and 

0)^^ aie found to interchan gp twice Y/ith each other when 0 = 45° v/hereas 
they do not interchange when 8=0®* As the fibre -orientation angle is 
inercasody the frequency curves for ^^2 ^21 th^se for 

cone closer and coincide for ^ = 45° • 

The variation of a> \7ith d/a is sliewn in Figures 3.15 and 3.16 
for a unidirectional reinforced boron-epoxy plate in the case of 8 = 0® 
and 8 = 45° respectively* ‘The frequencies increase 

with d/a ratios for both 0 = 0° and 8 = 45°. The behaviour of frequency ^ 

0)^2 iSj hoYrever^ different it increases with d/a ratios when 0=0° 

(figure 3.15) YThereas it is found to def^reaso with increasing values of 
d/a and then to increase at higher values of d/a Y/hen 8 = 45° (li-gure 3.16). 

The behavjour of frequencies ^21^ ^22^ ^14 ^23 

sinilar in case of fib re -orient at ion angles 0° and 45°. The frequency 
cuorvGS for °^21 those for ^23 closer with the 

increase in the fibre orientation and they coincide when 8 = 452 

For unidirectional graphite-epoxy plates 5 the variation of with 
d/a is shown in Figure 3.17 for = 0°. The frequencies ^^<1? ^15 

corresponding to symmetric inodes increase with d/a. The behaviour of 
frequencies ^14^ ^21 ^22 similar to that found for unidirec- 

tional glass--epoxy and unidirectional boron-epoxy plates. It is interesting 
to note that there is no in iBr change of the mo dos shown when 0 =45° 

(Figure 3.18). Here also^ as for 0 = 0° j increases with d/a. The 
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frequencies *^22 *^13 'fc'^ decrease with 

increase in d/a values^ the decrease being nore in ^^^2 coisparea to that 
in the fiase of G = 0’^ . Those frequencies, however, increase with 
higher values of d/a. 

b) The Effect of Pibre-Orientation Jingle on the Hataral Frequency 

The influence of fibre -orientation angle has been shov/n on fre- 
quencies and selected values of d/a ratios. (Figures 3.19 

to 3.21). For a. partic'’lar value of G ^ is found to increase with 

d/a ratios; tliis is not true for “^2^ frequency curves for d/a = 0,3 
always lie below the frequency curves for d/a = 0.0 and d/a = 0,5 in all 
the cases considered. 

c) The Effect o f Modulus Eatio on the Natural Frequency 

In order to understand the effect of modulus ratio on the natural 
frequency of plabcs with holes, the behaviour of for different 

modulus ratios is shown w?ith d/a in Figure 3.22 in the case of fibre 
orientation an^cs 0° and 45°. It is found that increases with 

d/a for an modulus ratios. The frequency curves for 0=0° and 0 = 45° 
for a given mociilus ratio, except in the case of balanced bidirectional 
plabes arc foixnd to intersect with each other. This behaviour was also 
observed in clampcd-clampcd plates f/ith cut-outs. 
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d) The Effect o f Hole on the tfo * Shape 

The eivonvoefcors normalized v/ith respect to the mass matrix: have 
been obtained in all the cases considered. Table 3.6 givos a comparison 
of the frequency values for a symmetric mode and the normalized eigen- 
vectors for all raodilus ratios for a given cut-out and hole parameters. 
The results are given for the fibre-orientation an^ e 6 = 0°. 
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CHAPTER 4 


EXPERKEl'JTAL STUDY OH VIBRATIOH OP PLATES T/ITH HOLES 


4.1 Introduction 

The method of analysis described in the previous chapters has been 
verified by cougjaring the results obtained v/ith available results for 
isotropic plates. It was felt desirable to compare the numerical results 
for ortho tropic plates with expci±mental values. Por this purpose, 
experimental work has been carried out to study the effect of holes on 
the natural frequencies and mode shapes of simply -supported and clamped- 
clamped square plates. Cut-outs are not considered in the experimental 
study because of the practicsil difficulties in making them. The experi- 
mental set-up and procedures have been described in this cbiapter. The 
experimental study consisted in characterizing the composite material, 
achieving the approximate boundary conditions (simple -support and clamped- 
clamped) by testing solid isotropic plates for which results are known 
and finally testing the composite plates. 

Humerical results have been computed for uni directionally 
rcinfbrced E-glass-epoxy plates by taking dynamic properties. The theore- 
tical results for a plate with hole are compared with experimental results 
and are also congsared with those for a plate with central square cut-outs. 
The theoretical mode shapes have been obtained and conpared with experi- 
mental results. The effect of hole size on the normalized eigenvectors 
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and iitac effect of anisotropy for a particular hole size on the natural 
frequency and corresponding eigenvectors arc also investigated in this 
chapter . 

4*2 Characterization of the Conposite r.-aterial 

In order to study the response of uni directionally reinforced 
composite structural elements under dynamic loa.ding conditions, it is 
necessary to study the dynamic behaviour of these materials. Various 
methods for determining dynamic stiffness and damping :of composite 
materials have been reviewed by Bert and Clary Dod- The methods that 
arc comncnly used for determining dynamic modulus and damping coefficients 
are 

1) firoe vibration technique, 

2) hysteresis measurements during sinusoidal excitation, 

3) forced vibration response technique, 

4 ) pulse propagation technique. 

Ihe last two methods are the most pop-ular methods in composite 
material dynamic -characterization. Itorced vibration response technique, I 
which is simpler, has been used here. In this method the mdulus is 
determined from the resonant frequency and damping from the measxirement 
of half-power point frequencies. The method consists essentially of 
exciting a thin cantilever beam (or any other specimen configuration) 
and using the peak amplitude as the criterion for determining the resonant 
frequencies. Thus the modulus and damping can be calctilated at each 
* A part of Chapter 4 is reported in [88, 89 j[ . 
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natural mode of vibration. The disadvantage of this method is tha.t when two 
frequencies are closely spaced, both the modes nay contribute to the ampli- 
tude 0 f an intermediate frequency resulting in a pseudo-node. 


For thin homogoneous isotropic beams with small internal damping 
(loss factor loss than 0.1), the Bernoulli-Suler tlieoiy of vibra.tions 
of clastic beams gives the relationship between the real part of the 
comploc dynamic modulus and the resonant frequencies of the beam as 


(V) 


(4.1) 


where g = gravitational constant 
h = depth of the beam 
Y = sp ecific weight of the beam 
1 = free loogth of the beam 

= corrected resonant frequency for the n-th mode 
= characteristic laode constant 

Bernoulli -Euler theoiy does not indudo the effects of shear 3 

deformation and rotatory inertia which arc important at higher modes of 
vibration. Ifereover, the shear deformation effects become significant oven 
at low modes of vibration in hi^ly anisotropic beams for which 
ratio is greater than 25 [jll]* present case E-ji /^12 is 3.®ss 

than 25 and therefore, only at higher modes the shear deformation effects 
become significant. 
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The corrected resonant freque-nqy, 


u , is given 'ey 
nc 


“ac 



ir^ , nr ^2 


(1 


_E_ 

KG 


-1 


V7here ^ 

n = tho mo do number 
r = radius of gyration 
G- = shear modulus 
E = tonsil G modulus 
K = shape factor constant 

= measured natuia.1 frequency 
^nc “ corrected natural frequency 

The damping ratio ? is given by 


(4.2) 


U2 ~ “i 

C = 

2 01 
n 


(4.3) 


where and oi^ are half-power point frequencies. 

The composite material used to study the d 5 ?namic response of 
composite plates was uni directionally reinforced E-glass epoay material, 
Moesurements of dyaamic motiilus and damping ratios were made on cantilever 
beans with fibre -orientations 0°, 15° j 30°, 45° and 90° with respect to 
the beam axis. The specimens wore cut from unidiroctionally reinforced 
sheets and machined to size on a high speed routing machine. The sheets 
were made by a hand 1^-up technique. The dimensions of the beams were 



1” widCj *094^' thick and 9*1*' long (The variation in v/idth VvT.s -f *001” 
and that in thickness was .0015”)* 

Figure 4.1 shows the vibration tost set-up which used rxi 13 niocki 
PM-50 exciter driven by an oscillator through a nodcl 2250 IB ariipliiier* 

The beams were mounted on an electro dynamic shaker# Sloctricril resistance 
foil gauges (33-4? AF-8) were fixed on top and bottom surfaces of the 
beams to record only bonding strains* This was done keeping in view the 
fact that tv/isting is not pronounced in ^ass-fibro reinforced eposy 
beams as reported by Schultz and Tsai Csd* 2’hc clamped end of the beams 
was given a sinusoidal harmonic motion and the output of the stra.in gauges 
was seen on an oscilloscope (Type 564B Teirbronix) with a carrier amplifier 
plug-in-*unit • The resonant; fi*equcncy was determined by seeing the peaking 
0 f the strain in the oscilloscope and then carefully adjusting the exciting 
frequency in order to get the maximum strain. The resonant frequency was 
measured on an electronic counter. The beams were excited over a 
frequenqy range of 10 - 6000 cps* The dynamic modulus was determined 
taking into account the effect of shear correction. A few cases were found 
where two closely spaced resonant peaks occurred instead of the norsml 
single peak response. These are believed to be associated with the 
stiffness of the damping arrangement, and v»?hen it occurred, the date for | 
that mode was discarded. 

The modal frequency ratios ^5^? 30*^, 45^ and 90® 

cantilever beams are given in Table 4.1 and are compared with the theore- 
tical value calculated on the assumption that the iBjdulus is frequency 
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independent and that the beam is a unifoiuly hoiaegGneous Bcmoulli-Enler 
beam. 

Figure 4.2 shows the variation of the c^ynanLc modulus with fre- 
quency for various fibre-orientations. The effect of shear correction 
term is shown only for a beam with fibre-orientation. The shear effect 
is seen to be pronounced only at higher mode frequencies. An average 
value of the modulus is chosen for nximerical computation. 

The sta.tic modulus was determined by testing tensile specimens 
of 0 ° , 45° and 90° fibre-orientations, with strain gauges mounted, on an 
Instron machine DU- 

Table 4.2 gives the dynamic (average) and static properties of 
\ini directionally reinforced E-glass-epoxy composite material used in the 
experimental study. The static moduli are seen to bo lower by 3-33 percent 
than the dynamic nsoduli. 

4.3 Vibration Study 
4.3.1 Test Plates 

This study was concerned with the vibration of simply -supported and 
clampcd-clamped plates with circular holes. The isotropic plates were 
made of aluminium, had a thickness of .048" and had dimensions of 
9.5" X 9.5". Eight isotropic plates were studied. The eight plates 
consisted of one solid plate and seven plates with hol^ of 0.45", 0.90", 
1.8", 2.7", 3.6", 4.5" and 5.4" diameter. Ihe test plates used in the 
experimental investigation other than alinainium plates were cut from a 
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Tatle 4.1 

llodal frequency ratios for a cantilever beam - 


I/ode 

Isotropic 

(iheoiy) 


G-ene rally 

orthotropic 

(measured) 


0° 

15^ 

30° 

45° 

90° 

1 

1 

1 

1 

1 

1 

1 

2 

6.27 

6.03 

5.91 

5.93 

5.95 

5.98 

3 

17.5 

17.11 

16.98 

16.98 

17.05 

16.76 

4 

34.4 


34.1 

35.2 

- 

- 

5 

56 .8 

55.0 

55.2 

57.2 

- 

54.5 

6 

84.9 

82.1 

81 .4 

82. 7 

85. 0 

82.5 

7 

119 

llt.t 

115.6 

117.6 

120.4 

117.9 

8 

158 

154.1 

154.2 

154.2 

165.1 

156.8 

9 

203 

195.2 



213.8 


10 

253 

248.4 





HZ 


23.23 

20.75 

16.80 

14.42 

14.13 


* (These frequencies (corresponding to blank spaces) could not be 
obtained because of the resonance of the clamping fixture or the 
presence of pseudo modes. 
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Sable 4»2 

A conparison of static and dynanic properties of a uni directionally 
reinforced E-glass-epoxy material , 

Static lynamic 

fi 6 

2.02 X 10 psi 2.08 x 10 psi 

0.53 X 10^ psi 0.74 x 10^ psi 

0.18 X 10^ psi 0.27 X 10^ psi 

0.35 

0.05 Ib/inch^ 
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uni directionally reinforced E-glass-epoxy sheet. Square plates with 
0° , 15°, 30° and 45° fibre orientations wore tested. The principal elastic 
constants of the conposite naterial wc-re found by forced vibration response 
method and has been discussed earlier in Section 4.2. In order to study 
the effect of holes on the vibration characteristics of uni directionally 
reinforced plates, the following steps were taken. 

a) Vibration clnracteristics for the unidirectionally reinforced solid 
plates were fcund . 

b) Small holes were drilled in the plates. 

c) The holes were enlarged to size on a high speed routing machine 
using aluminium templates. 

d) Vibration characteristics for the unidirectionally reinforced 
plates with holes were found. 

e) The holes were increased in size by machining again and the 
corresponding vibration chara.ct eristics for plates were found for 
each predetermined hole size. 

The unidirectionally reinforced plates were 9.5" x 9.5" square f 
having a thickness of .094". The hole sizes used in the test were 
0.45", 1.8" and 3.6". 

4.5.2 Simple Support and Clamping Arrangements 

A test fixture v/as designed and fabricated to stu(^ the vibrations 
of plates. It ck)nsisted of four vertical members bolted to Uie foundation 
and four horizontal members supported by the vertical members. On the 
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horizontal cenbors a thick plate with an internal square opening was kept. 

The thick plate could bo socurca on the horizontal menbers. Provision 
was nade on the thick plate to place rollers for conducting axpcrinents 
on simply -sup ported plates. Simple supported conditions were a.chioved 
by placing test plate bGtv;een two sets cf rollers (Figure 4.5). The upper 
set of rollers were held in position on another plate with an internal 
square opening. The test plate was held between rollers by bolts. Tests 
were carried out on a solid aluminium plate and the torque on the bolts 
to hold the thick plates with rollers was appropriately chosen to give a 
good agreement between the theoretical and experimental results. The 
same testing conditions were maintained throughout the course of experi- 
ment. 

For clamp Gd-clampcd conditions, the test plates were kept on the 
thick plate and clamped by means of two sets of thick mild steel flats by 
nuts and bolts. The effect of various torques applied for tightening the ) 
bolts mis seen on the natural frequencies of an aluminiuiii test plate without j 
hole. With increase in the value of torque applied on bolts, the frequency 
was found to increase. But after a certain value of the torque there was 
no app3recia.blG change in the frequency. A particular value of the torque ' 

(14 lb. ft) was chosen and maintained in testing all plates. Several 
tightening sequences did not produce any variation in the natural fre- 
quencies. The frequency values for the aluminium plate were checked with 
available results and a good agreemoit was found. This su^ested that the 
desired edge conditions have been obtained. The eseperinents on luiidirec- 
tionally reinforced E-glass-epoxy plates were conducted in identical conditions. 
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4 .3 #3 Source of Excitation 

The test plates were excited acoustically by cierns of a loud 
speaker (shov/n irx Figure 4.4)^ field of the coil was driven by a a. c. 

supply. Tiie signal was fed to the loud speaker fnorn an oscillator through an 
amplifier a.nd the amplitude of the signal could, be controlled both from 
the oscillator and the amplifier. The signal frequency was accurately 
measured by an electronic ccuntor. The test set-up is shown in Figure 4.5. 
Ihe loud speaker was used because of the following reasons: 

1) In the case of excitation by an elcctrodynaiiic shaker the whole 
fixture with test plate has to bo mounted on the shaker which may not 
be always possible duo to the force limitations of the shaker. Also^ the 
test fixture resonance may create difficulty in obtaining a nodal pattern 
when one of the natural frequencies of a test plate happens to be near 
the resonant frequency of the fixture. 

2) The variable frequency air jet excitation method can be used but 
the arrangement is quite complex as a continuous oir supply is always 
required. 

Ihe position of the speaker did not affect the natural frequencies 
of the plate. 

4.3.4 Measurement of Natural Frequency 

A contactless vibration pick-up (type m¥‘-200) was used to find the 
natural frequencies of the vibrating aluminium plates. Ihe pick-up 
generates an electrical output. The signal from the pick-up was 

Manufactured by the National Aeronautical Laboratory, Bangalore (India). 
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amplified and displayed on an oscilloscope. Ihe input signal was also 
fed on the oscilloscope and no change in the frequency v/as foxmd between 
the input and output signals. 

The first step in obtaining the natural frequencies, ?/hich in this 
experiment were indicated by the maxirium pick-up output, was to place the 
pick-up at a point above the plate. The frequency of the signal from the 
oscillator was slowly increased and the output of the pick-up observed 
on the oscilloscope. At a certain frequency the output of the pick-up is 
fOTjnd to be naxiiium. The frequency is carefully adjusted at this point 
so as to get the maximun output displayed on the oscilloscope and is then 
measured on an olcctronie counter. Nodal patterns were also used to measure 
frequencies and these agreed well with those given by the contactless 
pick-up. 

In the case of imi directionally reinforced plates, the natural fre- ^ 
quencios were measured when a clear modal pattern was observed on the 
vibrating plate. 

4.3.5 Nodal Pattern Detection 

When a plate is excited at a natural resonant frequency it vibrates 
in a particular mode shape. This shape can be visualized by sprinking dry- 
sand or aluminium oxide granules on the surface of the plate. The particles 
collect on the nodal lines where there is relatively less motion, "fliexi the 
plate vibrates at a frequency other than the rescaaant frequencies there are 
no nodal patterns and the granules bounce forming no pattern. With the 



Test plate 





a. Simple-support arrangement 
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b. Clamped -clamped arrangement 


Fig. 4-3 Simple-support and cfdmping arrangemeni 
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granules on the surface of iilie plate, the input signal frequency was slowly 
increased until the motion of the aJuniniun oxide granules indicated a 
resonant condition. The frequency v/as then tuned to produce the nost 
clearly dcfincxi nodal pattern. The frequency v/a.s reasured on the electro- 
nic counter and the node shape photographed. The upward frequency sweep 
vi?as then centinued until the next rescriance was obtained, the frequency 
tuned in and measured and the mode shape photographed. This was continued 
to get higher frequencies and corresponding lode shapes. 

Because of the loud speaker pov?er limitation only the first threc- 
symmetric EXJdes could be excited in the -unidiiectionally reinforced plates. 
However, these were considered sufficient for the verification of theore- 
tical results. 

4.4 Results and Discussion 

Natural frequencies and Doda shapes have boon obtained for square, 
simply-supported and clampod-clampod imidircctional glass-epoxy plates with 
central circular holes for different fibre-orientation angles. A comparison 
between the experimental a.nd theoretical frequencies for a sin 5 )ly-supportod 
and clampc3d-claE55Gd solid isotropic plate has been made to ensure that the 
required edge cenditions have been approximated. Itor a comparison of 
exp crimcntal and theoretical frequencies the computations for un-idirec- 
tionalHy reinforced 1-glass-epoxy plates with holes were made taking the 
dynamic properties from Table 4.2. Results are plotted in a non-dimensional 
form. 1 comparison of experimental and theoretical mode shapes has cilso been 
made for unidirectionally reinforced glass-epoxy plates. 
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4.4.1 Simulation of Edge Conditions 

Table 4.3 gives a comparison of the experimental and theoretical 
frequencies for a solid isotropic plate. The following properties have 
been used for the alumitiium plate tested. 

fi 

E = 10.5 X 10 psi 

pg = 0,096 Ib/inch^ 

V =0.3 

The agreement between the theoretical and experimental values is 
fairly good. The percentage variation in the frequency values for a 
simply -supported plate is maximum for the f-undamental frequency but this 
variation is found to decrease for higher modes. In the case of clamped,- 
clanped plates, the effect of clamping has been taken into account by 
considering the effective length as suggested by Laura C94II* They 
defined the effective length as, 

a^^^ = a + s/n' (4.4) 

where a is the exposed length of the plate, 

S is the clamped length, and 

n’ is a constant for a particular clamping. 

A proper value of n' was chosen so that the theoretical frequencies 
compare well with the experimental values for a solid isotropic plate. 

The value of n' , once fixed (2.34), was used for all the other experiments. 
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Table 4.3 

A comparison of theoretical and experimental frequencies (liz) 
for a solid isotropic plate . 


Eoundaiy 

Conditions 

Simply 

'-supported 

Claraped* 

-claiiiped 

i-fode 

fheorj 

ExperiiTient 

Theoiy 

Experiment 

1 

117 

126 

249 

240 

2 

291 

308 

508 

509 

3 

465 

482 

749 

- 

4 

582 

589 

911 

- 

5 

756 

772 

915 

92C 

6 

988 

- 

1142 

1150 

7 

1046 

1071 

1457 

1419 

8 

1510 

1476 

1526 

1580 
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It is vrorth Eientioning here that the effect of clcciping has not been 
considered by Ashton [jl^ also by ionahan {^S^U * reason 

that the theoretical frequencies obtained by then are higher than those 
obtained experinentally; the theoretical fundruiental frequency obtained 
by Ashton [jl a clamped-clanped boron-epoxy plate (0-90°) is approxi- 

mately thirty percent hi^er than the experimental value and the error 
in the fundamental frequency for a rectangular isotropic claEp)ed-clamped 
plate obtained by Monahan fourteen percent. 

4.4.2 The Effect of Discontinuity on the Uatural Frequency 

The effect of cut-outs and holes on the natural frequencies of plates 
has been discussed at length in Chapters 2 and 3. However, all the results 
have been shown for fibre -orientation nn^es 0° and 45° only. In this 
section, the variation of frequency with d/a ratio is shown for fibre- 
orientation angles 0°, 15°, 30° and 45° for uni directionally reinforced glass- 
epoxy plates. The theoretical results are compared with experimental results 
and are also compared with those for a plate having central square cut-outs . 

a) Simply -Supported Plates 

The variation of frequency with hole parameters d/a for unidirectional 
E-^ass-epoxy plates is shown in Figures 4.6 to 4.9 for fib re -orientation 
angp.es 0°, 15° » 30° and 45°; the solid lines show the variation of the 
frequency with hole parameter whereas the dotted lines show that with the 
cut-out parameter. The results have been shown for cut-out ratio upto 
0,8 whereas for plates vsith holes the frequency values coul d not be obtained 
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for some iiDdes beyond d/a = 0.4 due to computational diffieulties such as 
non convergence of the frequency values ■ and hence excess computation time 
involved) and large round-off erro3?s because of the numerical integration 
scheme uaad. 

The fundaiTiental fre. 5 .usnc 7 is found to decrease initinlly vjith d/c. 

ratios for fibre-orientations 0° and 15®» But as the fibre-orientation an^e 

is increased further, the rate of decrease of becomes less as shorn 

in Figures 4#8 and 4*9 for fibre -orientations 50® and 45^* The behaviour 

of frequencies m ^2 ^21 similar for all fibre -orientation angles* 

Small holes do not seem to influence these mcdes* However, with the increase 

in hole size, there is an abrupt decrease in the frequency values for these 

modes and they again increase as the hole size is increased. The frequencies 

Woo and w _ are found to decrease continuously with increase in hole size 
dc 25 

but these also increase at laiger d/a. ratios. For small hole sizes, the 
frequency is constant but then increases continuously as the hole size 
is increased for all fibre-orientations. There is a tendency for the modes 
corresponding to and to interchange with other for ail fibre- 

orientation angles whereas those corresponding to and interchange 
only in_case_qX,„J_ 50® and e = 45"^* As the fibre -orientation angle is 
increased the frequency curves for ^21^ ^25 ^32 closer 

and coincide for 6 = 45®* 

The experimental results corresponding to the first three symmetric 

modes (aT-.-,, and w" ) are compared with theoretical results for 

^ 11^ 15 51 

fibre-orientations 0®, 15® and 30®' and those corresponding to the first and 






0.2 0.3 0.4 0.5 O.S 0.7 0.8 

d/a 

d for 0 simply *" supported plate ; 8 = 15*^, R 1.0 



FlgA8 w vs <i/0 for a simply- supported plate;8 ~ 30 , R~l(^ 
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third syismefcric nodes for 8 = 45°. The agreement between the theoretical 
and experimental frequency values is good. Any descrepancy between the 
results can be attributed to the fact that it is extremely difficult to 
achieve perfect simple -sup port edge conditions experimentally. Also, there 
is a variation in the thickness of the plates (of the order of ^ 3 percent) 
which may give rise to a sligjit cha.nge in the modulus of the plates. 

A comparison of theoretical results between cut-outs and holes has 
also been shown in Figures 4,6 to 4.9. In general, it is seen that the 
variation of different frequencies for a hole and a cut-out are similrx. 

At smaller, d/a ratios (upto 0.15 - 0.2), the frequency values for a given 
mode in case of a cut-out as v;ell as for a hole are practically same. This 
suggests that for s m a l l d/a ratios, one can get the necessary information 
by computing frequencies for plates with cut-outs, the analysis for which 
is much simpler, 

b) Clamp ed-Clanped Plates 

The variation of frequency with d/a for unidirectional E-glass-epoxy 
plates is shawn in Figures 4.10 to 4.15 for fibre -orientations 0°, 15°, 30° 
and 45°. The theoretical results are cemparod with those for a plate with 
cenfral square cut-outs and also with experimental results. It is seen 
that for all fibre-orientations increases with increase in d/a ratios. 

This is not true for frequencies “■] 2 » *^21 ^23* constant 

for smaller values of d/a and decrease as d/a is increased. However, further 
increase in the hole size increases the frequencies. The frequency *22 
behaves in a different mmer; it decreases with increase in the hole size 
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but is foioiid to incxreasG with larger d/a ratios. For smaller values . 
of d/a there is no appreciable change in but it increases as the 

hole size is increased. For 8 =0'^ and 0 = 15° > nDdes corresponding 

to "tOgg “-13 sire found to interchange but these do not interchange 

for 8 = 30° and 8 = 45°. Ifodes corresponding to and inter- 

change only for 0 = 30° and 6 = 45°. It is interesting to note- that 
frequone;/ curves for ^^2 “qi closer and those for and 

*^22 from each other as the fibre -orLont at ion angle is increased. 

Experimental results for the first three symmetric modes are 
compaxed with theoretical results and the agreement is found to be good; 
here also, as in simply -seqsported plates, experimental results could be 
obtained only for the first and third symmetric nedo when 8 = 45° - A good ^ 
agreement between the cocperimental and theoretical results stress the 
need for considering the effect of clamping which is quite often ignored. 

In Figures 4.10 to 4.13, a comparLson of theoretical results 
between cut-outs and holes is sh>vm. IXae to oomputatiorGl difficulties, 
frequenqy values for son^ of the modes could not be obtained beyond d/a 
equal to 0.4 whereas results have been obtained for c/a ratios upto 0.7. 

It is seen from Pigu3res 4.10 - 4.13 that the vailation of different 
frequencies for a bole and for a cut-out are similar. At smaller d/a 
ratios, the frequency values for a cut-out as well as for a hole are 


practically sas®. 


13 O 


0.2 


m 


0.4 


0.5 


0.6 


d/a 

vs d/a for a clamped -clamped plate 
= 0®R = l.0 


m vs^ d/a for o clamped •‘ciarnped plate 
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4.4.3 IvSocle Shapes 

It becoinos easier to obtain the node shape for a particular 
frequency once the corresponding eigenvectors are known# The theoretical 
m d3 shapes have been obtained for both isotropic and iinidirectionally 
reinforced glass-epoxy plates. Only the node shapes for the conposite 
plates are shown and these are compared with experimental niDcle sho-pes. 

In obtaining the mode shapes for clamped-clamped edge conditions^ the 
computations have been carried out in double— precision arithmetic to 
avoid large round-off errors in function evaluation by making use of the 
values of coefficients given in Table 2*1 • The effect of the follo¥/ing 
parameters has been considered, on the node shapes of a utiidirectionally 
reinforced glass-epoxy plate for simple-support and clanped -clamped 
edge conditions. 

a . Fibre-orientation angle 5 

b# Presence of hole 5, 

c. Shape of the discontinniiy , 

d. Size of the hole. 

a. Fibre -Orientation Angle 

The theoreticp-l node shapes have been obtained for all fibre- 
orientation a.ngles when d/a — 0*0 and d/a = 0*2. However^ they are 
shown only for 0^ (Figures 4*14 ^d 4.16) and 45^ (Figures 4.15 and 4.17) 
fibre-orientations^ ♦ This is due to the fact that the node shapes for 
fibre-orientation angles 15® and 30^ are similar to those for fibre- 
orientation angle 0 ^ . The node shapes shown in Figures 4. 15 and 4*17 for 
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plefc es with fibre -orientation an^c- 45"^ haring simple-support and clamped 
edges respectively are similar to those found for corresponding isotropic 
plates. This has also been seen [^343 ^ simply -supported unidirec- 

tional boron-epoxy plate having fibre-orientation angle 45" where the 
authors obtained mode shapes which are in fact a linear combination of 
modes and are similar to these obtained in the case of an isotropic plate. 

A good treatment for the mode shapes resulting from a linear combination 
of modes has been given by Hikami |^5o3* experimental unde shapes t 

(shown for all fibre-orientation angles) are given only for plates with 
holes in Figures 4#18 and 4*19* A comparison of theoretical and experi- , 
mental mode shapes for fibre -oilen tat ion angles 0® and 45° for a. plate 
with a ix>le show good agreement. l}hG distortion in experimental mode 
shapes are due to l) the edge conditions being approximate and 2) the 
variation in the thickness of plates (the variation was of the order of 
+ 3 percent)* 

b. Presence of Hole 

The change in the mode shapes resulting from the presence of 
holes in plates is shown in Figures 4.14 to 4*19 for both simple -support 
and clamp ed-clamped edge conditions. It is observed that the symmetric 
modes get distorted due to the presence of bole in all the cases. This is 
due to the symnetric location of the hole with respect to the plate gcometr;}’^. 
Asymmetric location of the lx>le may effect the mode shapes for antisymmetric 
modes* The experimental mode shapes for plates with holes are shown in 
Figures 4.18 and 4.19 and they c 0 iiipa]i?e well with theoretical results. 
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Fig. 4. "18: Experimental Msde Shapes for a Simply-Supported Plate; 
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Por n siiaply -supported plate with fibre-orientation angle equal to 
45°, the presence of hole gives rise to a now mode; this has boon discussed 
earlier in Chapter 2. The interchange of fourth and fifth modes for a 
cl anpod-clampcd plate- can be seen in Figure 4.16 for fibre-orientation 
angl e 0° . 

c. Shape of the Discontinuity 

A coraparison of frequencies and mode shapes between a cut-out 
(c/a = 0.2) and a hole (d/a = 0,2) for simply-supported a.nd clamped-claiEped 
plates are shov/n respectively in Figures 4.20 and 4.21 for fibro-orientation 
angle 0°. fhe inde shapes aie found to be similar. This is true for other 
fibre-orienta.tion angles also. 

d. Size of the Hole 

It was found that the presence of holes distorted symmetric Esode 

shapes which show the coupling of different modes. For this reason, a | 

< 

coEipaiison of frequenqy values and mrmlized eigenvectors for a square, / 
sinply-supported and clanpod-clamped unidirectional plate is made for 
diffe3?ent d/a ratios for a particular mode. Table 4.4 gives this compari- 
son fbr fibre -orient at ion angle 0° and mode observed that as 

the iK>le size is increased, the contribution of other modes becomes 
more. 

Finally, the effect of modulus ratio on the ftequenqy values and 
normalized eigenvectors for a square, simply -supported and tlasped-elamped 
plate is shown in a tabular form (Table 4.5). The relative va.lues of 
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amplitude coefficients suggest that the- mo x shape is not iiuci: affectec 
by the modulus ratio for a given mode. 

4.4»4 Effect of Static and lynanLc I-S>duli on Eroquency 

The computations of natural frequencies have been done by making 
use of dynamic properties from Table 4.2, It v/ould be interesting to 
compare these frequency values v/ith those obtained using static piopearties. 
For this purpose, a comparison between frequency values computed using 
static and djmamic modxili has been made. The results are given in Table 4.6 
for a simply -supported imi directional glass-epoxy solid plate for three 
symmetric modes and the theoretical results are also compared with 
experimental results. It is observed that the theoretical values obtained 
using static nKDduli are lower than those obtained by using dynamic 
moduli. <"'hcn a comparison is made with experimental values, it is found 
that the error in the values of the first and third symmetric modes is 
comparativehy less than that for the second symmetric modoj the difference 
in theoretical (static mociili) and exp orimental values for the second 
mode is approximately sixteen percent ( 6 = 0°) whereas the difference 
in theoretical (dynamic moduli) and experimental values for the second 
laode is just one percent. 

The above observation led to the decision to use the dynamic 
properties in the computation of frequencies. 
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TRi'JJSISJT ?53?C::32 07 COISOSIIE PLATES 

5-1 Introduction 

I>ynamc analysis of the response of any structure or a structural 
element to shock and other transient forces is complex. Lany theoretical 
and numerical analyses of the ^mamic response of structures and structural 
elements have been made and these include classical normal mode solutions, 
lumped-parameter models, finite difference models, finite element models 
etc. These models are used to predict the dynamic response of structures 
subjected to transient loading. However, to develop a reliable dynamic 
analysis the following information is needed. 

1) the nature of the transient loading (e.g. peak pressure, or peak 
force, duration of the pulse, rise time etc.), 

2) the geometry and structural configuration of the element, 

3) the static and dynamic material characteristics of the structure, 

4 ) the (dynamic characteristics of the structure (e.g. natural 
frequencies) which is depen <tent on the boundary conditions, 

5 ) coinparison with experimental results. 

On the basis of the infomation needs enumerated above, 
theoretical and experimental studies were made on clamped- clamped 
rectangular isotropic and composite plates with and without holes. The 
theoretical stu<^ deals with the el^tic response of the plates uslr^ 
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classical rtorcal node approach and the convergence of the rcs’.ilts has 
been studied by taking different nunber of nodes. In the experinontal 
study, a shock tube has been used as the- ioadinc device-. The reproduci- 
bility of the pressure and strain neasurenients has been checked. !■>. 
conparison of the peak dynamic strains, experimental and theoretical, 
has bear nade for clamped-clanped plates. Ihc effect of holes on the 
dynazric response has also been studied by conparir.g peak strains for 
plates with and without holes. A good agreement is found for strain- 
time histories obtained experimentally and with theory. Finally, dynamic 
amplification factors (defined as ^ the ratio of the peak dynamic strain 
to the static strain) have been obtained and con^pared for isotropic and 
unidirectional glass-epoxy plates. 

5.2 Experimental Setnrp and Procedure 

Experimental work has been carried out to find the effect of 
shock loading on the response of isotropic and unidirectional glass-epoxy 
plates with and without holes. The experimental set-up (shown in Figure 5.1) 
and procedure have been described in this section. 

5.2.1 Shock Tube and Production of Shock Wave 

The stock tube is a long chaitoer of constant cross-section which 
has two chambers, the high pressure and the low pressure, separated by a 
diaphra^. The diaphra^ can be ruptured to produce a shock wave which 
passes down the tube. The shock tube used is of rectangular cross-section 
and is shown in Figure 5.2a. The higja pressure chamber is 6” x 4’* in 



cross-secticn ana 36" long, lbs low pressure- siSo has five soctions each 
28" long and having 6" x 4" cross-sc-ctions . The rfect''ngalar section is 
r*ade of four aluroiniur. plates. The different auctions on the lev; 
pressure side- are joined by neans o-f flanges and o::.dc- air tight with 
rubber gaskets provided in between tiaen. The high pressure and low 
pressure chambers are joined by -a circular flange with provision for two 
0-rings. On the high pressure side, an inlet valve for fillin' it with 
cOQpressed air and a Bourdon’s pressure gauge to read the pressure are 
pro vi (ted . 

The aligrcient of the shock tube is checked before starting the 
experiment. To produce a shockwave, the diaphragm (cellophane paper has 
been used hare) is kept between the two chaiibers and then both sides are 
joined. The high pressure charteer is pressurized by means of compressed 
air from a compressor (Pigujre 5.1). When the desired pressure as indicated 
by the pressure gauge is reached in the high pressure chamber, the 
valvera is clo^d and the di^hragm punctured at the centre by means of a 
hi^ resistance wire carrying high current. This causes the sudden rupture 
of the diaphragm and a shock vseve is produced. 

5.2.2 Eoasurement of Pressure 

It is necessary to know the peak pressure, duration of the blast 
wave, rise time etc., so that coaparlsoas can bo made between experimental 
and computed results. For this purpose, two pressure transducers 
(Kistler type 623 b/ 601A) - one raounted ori the shock tube near the low 
pressure end to trigger the oscilloscope and the other munted on a stiff 
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plate to know the exact pressure distribution on the teat plate were used. 


2bo stiff plate- was r.o’onted in tlic sane 
pressure traoosducer was roountel cn the s 
x'his was done for two roaserjs ; 


position as the- test plate and one 
tiff plate vdth pitiper aiaptois. 


a) 


the vibration of the plate- should net affect the signal fre.r. the 


transducer, 

b) the pressure inside the shock tube nay be different at different 

sections and therefore, it is necessaiy to know the pressure-tine 
history at the end section where the test plate is nounted. 


A typical pressure-tico histoiy is shown in Figure 5.3a. 


The pressure fron the transducer mounted on the shock tube v;as 
also recorded and it was found that the pressure distribution and the peak 
value of the pressure obtained from this transducer is the same as obtained 
from the one mounted on the stiff plate. As the least count of the Bourdon 
pressure gauge was 1 psi, it was felt desirable to record the pressure 
from the transducer mounted on the shock tube for each test to avoid errors 
resulting from reading the pressure on the Bourdai gauge. Bor a fixed 
pressure on the hi^ pressure clianber, hov/cver, the reproducibility of the 
shock wave was ^od. 

It is also necessary to check the nature of the shock wave- 
whether the shock wave is a plane stock wave or not at the end section of 
the low pressure side where the test plates were mounted. For this purpose, 
two pressure transducers were mounted on the stiff plate at arbitrary points 
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and the similarity of pressure-time records (Figure 5.3b) show that the 
shock wave reaching the end section of the low pressure side where the 
test plate is mounted is a plane shockwave. If the shock wave were not 
a plane wave, then it is ne*essaiy to knov/ its variation over the area of the 
test plate. As the pressure measurements can be made only at discrete 
points, it is likely that more errors are intro ciiceu in idealizing pressure 
distribution over the area of the plate and this may result in less 
accurate results obtained theoretically when the idealized pressure dis- 
tribution is used for coa^jutation purposes. 

The transducers used were of piezo -electrical crystal type which 
give a charge output proportional to pressure. The charge output is 
convetdjed to "voltage and amplified by means of a charge anplifier 
(Kistler-type 566 M3). The pressure signal from the charge aii 5 >lifier is 
fed into a storage oscilloscope (Tektronix -"type 564B)- The pressure 
measuring device has a rise time of a few micro seconds (3 ^micro -seconds) 
and a bandwidth of several mega-qycles. Since the rise time and duration 
of the pulse are of the order of 0,05 milliseconds and 7-8 milliseconds 
respectively, the transducers were sufficient for the purpose used. 

5.2.3 Ifeasurement of Strain 

The strain signals were found by means of electrical resistance 
foil g auges (.SB4. type PAP-12-12, S6) laounted on the test plates. The 
position of the strain gauges on plates with and without holes is shown 
in Figure 5*4. As the duration of the blast load is of the order of 
milliseconds, these gauges were quite sritable rir the measurement of 
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dynandc strains. 'The output fron the strain gauges were recorded on a 
storage oscilloscope using a carrier amplifier plug-in-unit (tjpe 
Tektronix); the oscilloscope being triggered by the signal from the 
pressure transducer mounted on the wall of the shock tube. Ht-l-t-brido'- 
circuit was used with a dunmy strain #iuge in the other am. The straiii 
gauges were calibrated with a standard calibration resistor (150 
provided in the carrier amplifier unit before taking any mcc-^urements , 

The strains and corresponding pressures were recorded for all the gauges. 
The strains were normalized with respect to the piessure to compaie the 
experimental results with computed results. 

5.2.4 Test Plates 

Aluminium and unidirectional reinforced E-^ ass-epoxy rectangular 
plates were used for blast response studies. Tests were done on plates 
with and without holes. Only one hole size was considered. Figure 5.4 
shows the test plates and also the location of the strain gauges mounted 

on them. 

5.2.5 Clamping Arrangement 

The plates were clamped at the end of the shock tube by means of 
two supporting francs. One supporting frame was rigidly fixed to the 
shock tube. The test plate was rigidly clamped to this supporting frame 
with the help of another frame and bolts. Equal tightening of the bolts 

was adbiaved by hbbxib of a 




'1 

i 



fpfMrH ; schematic sketch of the experimental set-up for shock studies 
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Fig. 5.3a: A Typical Pressure-Time History 



Pig, 5.3b: Pressure-Time Trace as obtained, at Two Arbitrary 
Points of the Plate 
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5.5 Method of Solution 

The equations of motion for composite plates with discontinuities 
were derived in Chapter 2 and studies were made on the free vibration charac 
terlstics of these plates by substituting F(x,y,t) = 0 in equation (2.14). 
In the presence of an external force, the equilibrium equations given in 

equation (2.14) axe 


R S ii 

1 I - 

ift=1 n=1 




^inn 


+ 


R 

I 

] 0 fc :1 


S 

I 

n=1 


(C^^ - ot^^) q, 

^ Ml 


= p 


m 


xon 


Pi 3 


(2.14) 


where 


ij id and P . . have been defined earlier in Chapter 2. 

“mn’ ^mn' °mn’ “mn ^Pij 

It is possible to express equation (2.14) in a matrix form as 


[mC 

{q} 

+ 


where, 

’ pi 

= 

mci 

- 

^mn 

V 

= 

cij 

m 

id 

~ “mn 

\ 

= 

q 

m 

f 


= 

^PiD ’ 

and P 


(i 

- 1) S + 

1 

= 

(m 

- 1 ) s + 


{P} I 


(5.1) 


(5.2) 


(5.3) 


S is the- maximum number of terms 


taken in y-direction. 
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Mathematically, equation (5.1 ) represents a system of linear 
differential equations of second order and, in. principle, the solution to 
these equations can be obtained by standard procedures for the solution 
of differential equations with constant coefficients p 95 ? 96”! • However, 
the procedure proposed for the solution of general system of differential 
equations becomes cumbersome if the order of the matrices is quite large 
and the equations are coupled. The procedures that are commonly used to 
solve equation (5.1) are divided into two methods of solutions direct 
integration and classical normal mode approach. These two methods of 

solution have been discussed in detail in p97_|. The normal mode approach 

h 

has been used in the present study and the effect of damping has been ). 
neglected. 

5.3.1 normal Mode Method 

The mode superposition method consists in the transformation of 
djrnamic equilibrium equations into a more effective form using the follow- 
ing transformation on the generalized coordinates q^ 

’ ^ 5 . 4 ) 

where ]~G~j is a square matrix and (n) is a time -dependent vector of the 
same order as that of {q } ; the order of { q } is assumed to be n. The 
transformation is unknown and is obtained using the displacement solutions 
of the free vibration equilibrium equations with damping neglected (even 
if damping is present), 

pipi {q } + CKli{q} «= 


{ 0 } 


(5.5) 
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The solution of eq.uation (5.5) can be assumed as. 


{q.} = {A } e 


(5.6) 


where{A } is a vector of order n. Substituting for{ 1 } in equation (5.5), 


the general eigenvalue problem is obtained. 




'LmD tA> 


(5.7) 


solution of the eigenvalue problem in equation (5.7) yields 


2 2 

eigenvalues > (<^2 ■ 


and corresponding eigenvectors 


A }25 ... (a) . These eigenvectors are orthonormalized with 

respect to the mass matrio^jFj ; 

Clf h = 1 ; ' i = 3, 

= 0 ; i / D» 

2 2 2 2 
and 0 1 (0^ 1 ^^2 - ‘^3 *”• -“n 


(5.8) 


(5*9) 


The matrix [^(^3 defined as the mo.dal matrix whose columns are 
the eigenvectors {A}^. 7fe also define a matrix which stores the 

eigenvalues on its diagonal. Thus, 


^G-^] = Q{A},^j{A}2> ...j ^A}^ 


CO 


(5.10) 


2 
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Uow, the n solutions of equation (5.5) can be written as 

noil'll] = 

Premultiplying equation (5.11) by the transpose of jj^ H > S®"*" 

[E'railD = DO" Oil DO D^ 

Making use of the orthogonality relation in equation (5.8), 


we have 


nf ra HD - c« 3 

The transformation on the generalized co-ordinates given in 
equation (5.4) leads to the following equilibrium equations in modal 
gen oralized co-ordinates j 




(5.14) 


where , 


(5.15) 


The initial conditions (i.e. at time t - O) on 


' {Ti(t) } are given as 


{n(0) } = CC1"D®^0 ^ ’ 

«(o)} = i:([f'Dii(i(o)i . '' 

Equation (5.14) represents a system of n uncoupled differential 
eq^uations in modal generalized co-ordinates, {h (t)} with the initial 
conditions in equation (5.16) and can be written as 
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with initial conditions as 

Ti^(o) = {r}^CM3^9(o)> » 

^ 1 ^( 0 ) = ^ • ( 5 . 17 ) 

Applying Laplace transform technique, the solution of equation (5.17) 
is given as 

sin (o) t) 

n^(t) = n^(o) cos (u^t) + n^(o) — — 

t 

+ — / H (t) sin u) (t - t) dT ; 

u) P, r ^ 

r 0 

r = 1,2,..,, n . (5.18) 

If the modal generalized force W^(t) is known, equation (5.18) 
can be solved and p^(t) evaluated. 

5.3.2 Evaluation of E^(t) 

The method outlined in Section 5.3.1 is quite general and can be 
applied to any structural system. In the present section, the expressions 
for nodal generalized forces are developed for the idealized experimental 
pulse shape shown in Eigure 5.5a. 

The modal generalized force givsn in equation (5.15) is 

{E(t)} = 
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and therefore. 




(5.19) 


Wow, 

P = I’t,- • 
p PiD 

and from equation (2,13) > 

P(x,y,t) <t)^^ dA (5.20) 

It was discussed in Section 5.2 that the shock ware reaching the 
plate is a plane shock wave and is therefore, uniformly distributed over 
the surface of the plate. The distributed load is assumed to be of unity 
value and its variation with time is given as 


P(-fc) 


0 
1 

1 - 


t - t^ 
tg- I:-] 


_a> < t < 0 

t < t^ 


1^-1 


and t > tr 


(5.21) 


Thus W (t) can be evaluated from equations (5.2l), (5.20) and (5.19). 


5.3.3 Determination of Displacements and Strains 


In the determination of displacements and strains, it is necessary 
to develop the escpressions for generalized co-ordinates, S' 

coordinates q^ and the modal generalized coordinates are related and 

given in equation (5.4) as 
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and \ = V 

It is possible to evaluate ri^(t) from equation (5.18) once the modal 
generalized forces are known. The method for evaluating 1^(1:) has been 
given in Section 5.3.2. Substitutioii of equation (5.4) gives 

the generalized coordinates. The displacements and strains can be then 
found from the following expressions 


w(x,y,t) = 


e3j(x,y,t) 

ey(x,y,t) 


I I 

m=1 n;=1 


II 

1 

E . 

I 

nt=1 

S 

I 

n=1 


2 

9 4 (x,y) 

mn 

2 

9x 

I 

II 

E 

1 

n^l 

S 

i, 


2 

a 4 i (x,y) 
^mn 

2 


( 2 . 6 ) 

( 5 . 22 ) 

(5.23) 


5.4 Eesults and Discussion 

The method of solution given in Section 5.3 to find the transient 
response of a composite structure is quite general and can be applied to 
ar^r given structure if the admissible functions are piwperly defined for 
the structure . The numerical computations have been done for isotropic 
and con^osite plates in the case of sinply -supported and clamp ed -clamped 
boundaiy conditions. As experiments were conducted on clamped-clamped 
rectangular isotropic and unidirectional ly reinforced glass-eposy plates 
(6 = 0°), these results have been compared with the numerical results 
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obtained for claurped- clamped rectangular isotropic and unidirectional 
glass-epoxy plates; the material properties for the glass-epoxy plate 
being taken from Table 4.2. The classical normal mode approach was 
preferred because the computer program was already developed to compute 
natural frequencies and normalized eigenvectors. The reliability of the 
theoretical results computed depends to a great extent on the ^ccuracy in 
the measurement of the transient loading* This fact was observed by 
Baker et.al. Qdl] when they compared the theoretical elastic response of 
cantilever beams predicted from the linear elastic beam theory with measured 
blast response and emphasized the need for the accurate measurement of the 
transient loading* This fact has been kept in mind and therefore, the 
pressure-time history was recorded in each of the strain measurements. The 
theoretical strains have been computed for a unit value of the idealized 
blast load., shown in Figure 5.5a* The effect of the contribution of 
different modes as the dynamic response is studied and given in Table 5.1. 
for an isotropic clamped -clamped solid plate. In order to compare the 
theoretical and expeitlmental strains^ the latter have been normalized 
with respect to the pressure. The dynamic strains are also con 5 )ared with 
static strains. The theoretical and experimental strain-time histories are 
shown for a particular gauge in the case of a unidirectional glass-epoxy 
solid plate. Finally, the influence of holes is seen on the dynamic 
amplification factor for an isotropic as well as a composite plate* 
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5.4.1 The Effect of the Number of Modes on the I^smamic Response 

It is necessaiy to know the number of modes to be considered to 
get sufficiently accurate results for the ^mamic response of a structure. 
Eor this purpose 5 deflections and strains were computed for an isotropic 
and a unidirectional glass -epoxy clamped-clamped plate taking different 
number of modes and the results are given for central deflection and 
strains (1H, IV, 2H, ... etc.) in the case of a rectangular isotropic solid 
plate in Table 5.1. It is observed that a good estimate of the central 
deflection and strains 1H, 2H, etc. (strains in the x-direction) can be 
obtained just by considering only one mode- the fundamental mode. However, 
this is not true for strains IV, 2V, etc. (strains in the y-direction) . ^ 

Obtaining the dynamic response by just taking the contribution of only 
one mode may lead to absurd results. But as the number of modes taken is 
increased, good convergence is obtained and therefore, the results were 
computed by taking nine modes±H in all the cases considered. The same 
observation is found to be true for unidirectional glass -epojy plates. 

The contribution of hi^er modes becoming significant was also brought to 
light by Crocker ^ simplified analysis of the multimode response 

of panels, subjected to sonic booms. He mentioned that as b/a ratio is 
decreased (b < a) , the hi^er mode contributions to strains in the 
x-direction (i.e. strain in direction parallel to the longer dimension of 
the plate) become larger. In the present case, a < b and therefore, it is 
expected that the contribution of higher modes will be impoitant to strains 
in the y-direction and this has been found true. In the case of a square, ^ 
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isotropic plate axone expects that the higher mode contributions will not 
be significant as b = a ; however the same is not true for a square 
composite plate because of the orthotropy present. 

5.4.2 Comparison of Theoretical and Experimental lynamic Strains 

A comparison of theoretical and experimental dynamic strains is 
essential to ensure the reliability and validity of the theoretical 
model developed J^q^for the dynamic analysis of a structure. A comparison 
of theoretical and normalized experimental peak dynamic strains has been 
made for an isotropic and a unidirectional glass-epoxy composite plate 
and the results are given in Tables 5.2 and 5.3. The theoretical and 
experimental results for solid plates are also compared with those for 
plates having holes. The static strains have been given in all the cases. 

It is found from Tables 5.2 and 5.3 that the agreement between the 
theoretical and experimental peak strains for solid isotropic and unidirectional 
glass -epoxy plates is good except in the strain 3H for a glass-epoxy solid 
plate where the gauge did not seem to respond properly. The same is not 
true for plates with holes. The agreement between the theoretical and 
experimental peak dynamic strains is not very good; the maximum difference 
being 100 percent for strain 2V in the case of an isotropic plate. This 
suggests that the theoretical anabysis developed based on the classical 
normal mode approach is adequate for solid clamped— clamped plates whereas 
there is a need for more investigation into the transient response 
analysis of plates with holes. 
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The theoretical and experimental strain-time histories are 
compared in Figure 5.5b Sdt the strain gauge 2H in the case of a unidirec- 
tionally reinfoTced glass-epoxy plate. The agreement is good. The experi- 
mental peak strain^ however j is less than the theoretical peak strain and 
there is a time lag betv/een the experimental and theoretical strains after 
the application of the blast load. This may be due to the damping in the 
composite plates and also in the supports, 

5^4.3 Comparison of Dynamic Amplification Factors 

The theoretical dynamic amplification factors (defined as the 
ratio of the peak dynamic strain to the static strain) have been given in 
Table 5.4 for an isotropic and a -unidirectional glass-epoxy plate. The 
dynamic amplification factors for a solid isotropic plate lie in between 
the values 2.0 to 2,5 whereas those for a solid glass-epoxy plate lie in 
between 1.3 to 1,8 for the strains considered* This indicates that less | 
damage will be done to a glass-epoxy composite plate than an isotropic j 
plate when both are subjected to similar transient loads. The dynanlc 
amplification factors for plates with holes are approximately the same ' 

as the values for the solid plates, both in the case of isotropic and 
composite materials except in the case of strains 2"^ where the amplifi- 
cation factor is voiy high (6 to 7) w M(A seems unusual. An explanation, 
however, can be given if one compares the static strains for plates with 
holes and without holes. The values of static strain for 2Y is consi- 
derably reduced in the presence of a hole; this is justified as 27 is 
normal to the hole and nearer the hole geometiy. The' (^smamLc strains, 
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on the other hand, are not mch affected and therefore high values of 
dynamic amplification factors are obtained » 
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CHAPTER 6 


conclusions 


6.1 CondLusions 

The following are the major conclusions based on the present 

work. 

1) A method of analysis to study the free vibration characteristics 
and the transient response of composite plates with cut-outs (or holes) 
assuming the effect of cut-out (or hole) to be equivalent to a displace- 
ment-dependent loading has been developed and is found to give more 
accurate results compared to the other methods available. 

2) The frequencies for each mode in the case of a simply-supported 
plate with a cut-out (or a hole) vary with cut-out parameter (or hole 
parameter) in different ways for all the modulus ratios and do not uniformly 
increase or decrease, Por clamp ed-clamped plates also, the frequencies for 
each mode except the first vary with c/a (or d/a) in different ways, 

3) The fundamental frequency for simply -supported orthotropic plates 

(fibre-orientation an^e 9= 0°) decreases for medium cut-out (or hole) 

parameters unlike the behaviour in the case of isotropic plates wheare the 
fundamental frequency d e creases for small cut-out (or hole) parameters. 

The rate of decrease with c/a (or d/a) increases with the increase in the 



209 


modulus ratio. In the case of clamped-clamped plates, the fundamental 
frequency increases with the cut-out (or hole) size for all modulus ratios 
and all fibre-orientation angles; this behaviour is also fotmd for simply- 
supported plates when the fib re -orientation angle is 45°. 

4 ) The fundamental frequency for simply -supported solid plates 

(c/a or d/a ^ O.O) increases vdth increase in the fib re -orientation angle. 
The fundamental frequency, however, is found to behave in an opposite 
manner for a claiped-clamped solid plate; it deci-eases with increase in 
the fibre-orientation angle. This behaviour is observed for all modulus 
ratios censidored. Thus, increase in the fi.bre-o rientation an^e make ' 
simply-supported plates more stiff whereas it decreases the stiffness 

I 

of a claup ed-clamped plate when there is no cut-out (or hole). 

5 ) The behaviour of frequencies for balanced bidirectionally reinforced 
plates and for an isotropic plate is found similar. The fundamental fre- 
quency for a balanced bidirectional glass-epo^qy, simply -supported plate 

is lower than the fundamental frequency for an isotropic plate when 6=0° 
and is found to be higher when 6 = 45°. Tbi’ clamped-clamped boundary 
conditions, the fundamental frequency for an isotropic plate is found to be 
always higher than that for a balanced bidirectional composite plate for all 
cut-out (or hole) parameter considered. 

6) A similarity is observed between the behaviour of x frequencies 

b etween isotropic and the composite plates with fibre-orientation an^e 
equal to 45°. Tbr an isptropic plate, = “ 2 I ’ “23 ~ ~ ^^41 

etc., (frequencies corresponding to a combination of symmetric and 
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antisjnmnetric modes). In composite plates^ as the fibre -orientation 
angle is increased, frequency curves corresponding to ^^2 ^21^ 

and oT etc. a come closer together and coincide when 9 = 45®; this 

23 52 

means 

7 ) Two interesting features have been observed due to the presence 
of a cut-out (or hole) in a plate: a) the existence of a new mode in 

the case of a simply -supported plate and b) the cross-over of frequencies. 
In the case of a simply— supported plate, it is observed that the fourth 
mode is identical with the £Lfth mode for a solid isotropic plate* This 
is also observed for frequencies corresponding to both symmetric and 
antisymmetric modes. The introduction of a small discontinuity seems 
to perturb the system for which an exact solution is known and gives 
rise to a new mode. In the case of a clamped- clamped isotropic plate a 
closed form solution is not known and the approximation in the solution 
can be assumed to give rise to the new mode- for example, fourth mode - 

i 

even when there is no mit-out (or hole) in the plate. It is found that 
there is a tendency for the hi^er modes to interchange when there is a 
cut-out (or hole) in the plate. Mode shapes obtained at the cross-over | 
and nearby points for the frequencies which cross-over shov/ that at the 
cross-over frequencies the plate can have the mode shape corresponding to 
either of the modes which cross over. 

8) The presence of holes or cut-outs influences the modes in general. 
The influence on the modes depends on the location of the cut-out (or hole) 
in the plate 1 symmetric location of the cut-out (or hole) is found to 
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influence the ^mmetric modes and it is expected that asymmetric location 
will InfluGncG the antisymmetric modes. The presence of holes ^ therefore, 
distorts the mode shapes and this has been observed when experiments 
were conducted an isotropic and unidirectional glass-epoxy, square plates 
with centrally located circidlar holes 5 the symmebric modes were influenced 
most as the holes were centrally located. 

9) In general, the behaviour of different frequencies for a cut-out and 
for a hole is similar 1 The effect of the shape is seen only at larger 
cut-out ( or hole) ratios. 

10) An exparimental study was done on the free vibration characteristics 

of aluminium and unidirectional glass-epoxy, square plates. The unidirectional 
glass-epoxy plates v^cre made by a hand lay-up technique and the static and 
dynamic properties for the composite material was determined. The dynanic 
material properties were used in the numerical computation. A good 
agreement between the theoretical and experimental results has been 
obtained for isotropic as well as for the unidirectional glass-epoxj" 
plates. This proves the adequacy of the present method of analysis and 
stresses the need a) to use dynamic properties in the numerical computation 
of theoretical results to get reliable results and b) to take into account 
the effect of clamping which is commonly ignored. 

11) The mode shape for a particular frequency (identified by the largest 
amplitude coefficient in the corresponding eigenvector) is found to be 
independent of the modulus ratio for simply -supported and clamp ed-clamped 
plates. 
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12) An experimental study was mde on the blast respoiiso of isotropic 
and unidirectional glass-epoi^y rcctangularj cla up ed~ clamped plates using 
a shock tube as the loading device. The experimental results for strains 
have been compared with the theoretical results obtained using classical 
norroal mode approach for the isotropic and unidirectional glass-epo^^y 
c lamp ed-cl amp ed plates mth a.nd without holes* A good agreement' between 
the experimental and theoretical results for plates without holes show 
that the classical normal mode approach to the solution of transient 
elastic response problems is satisfactor;^^ for plates without holes* 

6*2 Recommendations for Puture \/ork 

Considering the direct extensions and generalizations of the 
present study pertaining to plates ^ the following problems are worth 
studying. 

1) An extension of the present analysis can be made to include the 
effect of bending and stretching coupling in the composite plates. 

2) Ihe present study can be generalized to include the effect of 
shear deformation and rotary inertia in composite plates with disconti- 
nuities . 

3) The effect of the cut-outs or holes can be studied on the dynamic 
response of composite shells. 
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4) The effect of damping which has been neglected in the present 
investigation can bo included in the free and forced vibration response 
analysis of composite plates and shells with discontinuities, 

5 ) Nonlinear vibration and nonlinear shock response studies can be 
32 B.dc on stnictujBl elGmcnts v/itli discontinuities. 
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INTEGRATION CONSTANTS FOR CDAIffi’ED-CDAMPED PIATES 



Defining X = x/a, Y = y/b, we have, 
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